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Abstract 

Massive black holes (BHs) are at once exotic and yet ubiquitous, residing in the centers of massive galaxies in 
the local Universe. Recent years have seen remarkable advances in our understanding of how these BHs form and 
grow over cosmic time, during which they are revealed as active galactic nuclei (AGN). However, despite decades of 
research, we still lack a coherent picture of the physical drivers of BH growth, the connection between the growth of 
BHs and their host galaxies, the role of large-scale environment on the fueling of BHs, and the impact of BH-driven 
outflows on the growth of galaxies. In this paper we review our progress in addressing these key issues, motivated by 
the science presented at the "What Drives the Growth of Black Holes?" workshop held at Durham on 26*-29* July 
2010, and discuss how these questions may be tackled with current and future facilities. 
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1. Introduction 

One of the most astounding astronomical discoveries 
of the last ^ 10-20 yrs is the finding that massive galax- 
ies in the local Universe host a central massive black 
hole (BH; Mbh ~ 10^-10"' Mq) with a mass propor- 



tional to that of the galaxy spheroid (e.g., Kormendy & 



Richstone 1995 



Magorrian et al. 1998 Ferrarese & 



Merrittl 120001 IGebhardt et al.! '2000' 'Tremaine et al.! 



2002| |Marconi & Hunt, 2003, Gultekin et al. 2009 )Q 
The tightness of this BH-spheroid mass relationship 
suggests a symbiotic connection between the formation 
and growth of BHs and galaxy spheroids. Identifying 
the physical drivers behind the BH-spheroid mass rela- 
tionship is one of the major challenges in extragalactic 
astrophysics and cosmology. 

BHs primarily grow through mass accretion during 
which the central source is revealed as an active galactic 
nucleus (AGN; e.g.,. Salpeter 1964, Lynden-Bell„1969l 
Shakura & Sunyaev|1973||Soltan|1982[|Rees|1984| i7A 
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' The evidence for this BH-spheroid mass relationship comes from 
tight connection between Mbh and the velocity dispersion, luminos- 
ity, and mass of the galaxy spheroid for x 50-100 galaxies in the local 
Universe. These constraints imply Mbh ~ (0.001-0.002) Mspi,. 



huge amount of energy is liberated during these mass- 
accretion events (e ^ 0.05-0.42 of the rest mass is con- 
verted into energy, depending on the spin of the BH; 
e.g., [Kot1[T963, Shapiro & Teukolsky|[T983] l, yield- 



ing large luminosities from comparatively modest mass- 
accretion rates: 



lO^^ergs- 



(1) 



where e is the mass-energy efficiency conversion (typi- 
cally estimated to be e » 0.1; e.g., Marconi et al.|2004 



Merloni|2004 1 and Lboi is the bolometric luminosity of 
the AGN. For example, the accretion of just ^ 1 yr"' 
(equivalent to the mass of the Moon each second) is suf- 
ficient for an AGN to outshine the entire host galaxy. A 
broad range of AGN activity is found even in the lo- 
cal Universe, with mass accretion rates ranging from 
^ 10^^-1 Mo yr-'. In Fig. 1 we illustrate some of the 
diversity in the AGN population within just ^ 300 Mpc; 
we note here that Sgr A* is not strictly an AGN but has 
been included to demonstrate that even BHs in relatively 
quiescent galaxies are growing to some extent. 

A direct connection between the growth of the BH 
and galaxy spheroid might be expected since both pro- 
cesses are predominantly driven by a cold-gas supply, 
which is provided by the host galaxy or the larger-scale 
extragalactic environment. However, the nine orders of 
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Figure 1 : The Eddington-ratio versus black-hole mass plane illustrating part of the range of mass-accretion phenomena seen in the local Universe 
{d < 300 Mpc). The plotted systems include radio-loud AGNs (Centaurus A, Cygnus A, and M 87), an ultra-luminous IR galaxy (Mrk 231), as 
well as typical low-to-moderate luminous AGNs (NGC 1068, NGC 4051, NGC 4151, NGC4258, and NGC 4945). All of the systems plotted 
here have reliable black-hole masses that are measured independently of the host-galaxy properties. The Eddington ratio of Sgr A, the accreting 
black-hole at the centre of the Galaxy, is « 10"*-10"' and cannot be directly plotte d on this figure. Th e curves indicate the space density of 
BHs growing at a given Eddington ratio, based on the synthetic AGN model of Merloni & Heinz] j2008} : <D = 10"' Mpc"^ (solid curve) and 
<I) = 10"' Mpc"^^ (dotted curve). The dashed line indicates a mass-accretion rate of 1 Mq yr"', which con'esponds to a bolometric luminosity of 
Lbol ~ 7 X lO*"* erg s"', which is equivalent to Z.2-iokcV ~ 3 X lO"*^ erg s"'. Systems lying below the solid curve are common mass-accretion 
events in the local Universe while systems lying above the dashed curve are rare events in the local Universe. Black hole mass and luminosity 
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[Braito et a l. 2004 , Y oung et al.| [2002 Di Ma tteo et al. |2003^ . ImageTredits: NGC 4945-2P2 Team, WEI, MPG/ESO 2.2-m Telescope, La Silla, 
ESO; NGC 706S-Francois and Shelley Pelletier/Adam Block/NOAO/AURA/NSE; Mrk 2J/-NASA, ESA, the Hubble Heritage (STScI/AURA)- 
ESA/Hubble Collaboration, and A. Evans (University of Virginia, Charlottesville/NRAO/Stony Brook University); Cygnus A-NASA, ESA, Neal 
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difference in pfiysical size scale between the BH and 
the galaxy spheroid would appear to preclude a direct 
causal connection. Many processes have been pro- 
posed which could forge a direct connection between 
the growth of the BH and galaxy spheroid, including 



galaxy major mergers, star-formation winds, and AGN- 
driven outflows. There is clear observational evidence 
for many of these processes occurring in individual sys- 
tems. However, it is often far from clear how univer- 
sal they are and what impact they have on the over- 
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all BH-spheroid growth. There are also many differ- 
ent observational and theoretical studies providing ap- 
parently contradictory results on the physical drivers of 
BH growth, with theorists often disagreeing with ob- 
servers (and observers/theorists disagreeing with other 
observers/theorists). A cause of these discrepancies 
may be that (apparently) contradicting studies are ex- 
ploring systems with different ranges in BH mass, Ed- 
dington ratio, redshift, or environment. 

To address these potential conflicts we organised the 
"What Drives the Growth of Black Holes?" workshop 
at Durham on 26*-29* July 2010. This workshop ex- 
plored the processes that drive accretion onto BHs, from 
the most luminous distant quasars to more quiescent lo- 
cal systems. A key aim of the workshop was to clar- 
ify the ranges of parameter space that are probed by 
different studies to better help understand how various 
key physical processes vary with these parameters. The 
workshop was organised into four main sessions that ad- 
dressed the following key issues: 

• How does the gas accrete onto black holes, from 
kilo-parsec to sub-parsec scales? 

• What are the links between black-hole growth and 
their host galaxies and large-scale environments? 

• What fuels the rapid growth of the most massive 
(and also the first) black holes? 

• What is the detailed nature of AGN feedback and 
its effects on black-hole fuelling and star forma- 
tion? 

One hundred and twenty participants attended across 
the 3.5 days duration of the workshop, the vast ma- 
jority of which presented scientific results: 53 gave 
oral presentations [18 invited; 35 contributed] and 49 
gave poster presentations j^j In this paper we review our 
progress in addressing these key issues, motivated by 
the science presented at the workshop. We also provide 
background material on the challenges faced in address- 
ing these issues to motivate discussion. However, given 
the large breadth of science covered by this review, we 
cannot provide a complete assessment of each paper in 
the literature for every component of every issue (see 
'Brandt & Hasinger '2005 ; Veilleux et al. ''2005'i Remil-" 
lard & McClintock 2006; Done et al..,2007, McNa-, 



mara & Nulsen 2007 , Ho 2008 , Shankar|2009[|Volonteri| 
|2010[|Brandt & Alexander|2010| for some recent indepth 



reviews on individual components within individual is- 
sues). Our aim is therefore to provide a solid ground- 
ing in the current scientific picture and hence a starting 
point for more detailed future investigations. 

2. How does the gas accrete onto black holes, from 
kilo-parsec to sub-parsec scales? 

The growth of the BH relies on the accretion of cool 
gas originally on scales orders of magnitude larger than 
the BH accretion disc, either from the host galaxy or the 
extragalactic environment (e.g., a companion galaxy; 
the dark-matter halo; cluster gas). The goal is to de- 
termine what physical processes deliver the gas from 
a: 10 kpc host-galaxy scales down to the BH accretion 
disc at r < 0.1 pc - an epic journey by any stretch of the 
imagination! Substantial barriers against the gas reach- 
ing the central regions are angular momentum and the 
gas collapsing and forming stars rather than accreting 
onto the BH. In this section we explore the mechanisms 
that can deliver gas into the nuclear region, including the 
competition for gas between star formation and AGN 
activity, the connection between accreting gas/dust and 
obscuration, and the process of mass accretion onto the 
BH. 

2.1. Driving the gas into the vicinity of the black hole 

The vast difference in size scale between the host 
galaxy and the BH means that the gas has to be driven 
down to a: 10 pc before it will come under the gravi- 
tational influence of the BH. The formidable force that 
needs to be overcome to deliver the gas from the host 
galaxy into the vicinity of the BH is angular momentum: 
the gas has to lose a 99.9% of its angular momentum to 
go from a stable orbit at r = 10 kpc down to r = 10 pc 



(e.g., Jogee 2006[ ). The challenges are to map the gas 
inflow from the host galaxy into the vicinity of the BH 
and to identify the triggering mechanisms of AGN ac- 
tivity. 

Large-scale gravitational torques, such as those pro- 
duced by galaxy bars, galaxy interactions, and galaxy 
major mergers have the potential to remove significant 
amounts of angular momentum and drive the gas into 



the central regions of galaxies (e.g., Shlosman et al. 



11989 1990[ Barnes & Hernquist 1992 


1996 Mihos & 


lHernquist"1996VBournaud & Combes 


2002| Garcia- 



The presentations can be found at the workshop web page: 
|http://as tro.dur.ac.uk/growthofblackholes/index.php 



|Burillo et al.^ 2005 l; see Fig. 2. However, the gravi- 
tational torques exerted on the gas by these large-scale 
processes will have a limited effect on sub-kpc scales, 
and smaller-scale processes (e.g., nested bars and nu- 
clear spirals; Englmaier & Shlosman|[2004 Maciejew 
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be found in all late-type systems (e.g.,|Laine et al. 



merger 



secular 



hot halo 



Figure 2: Schematic diagrams to illustrate the large-scale processes 
that are thought to be responsible for triggering AGN activity: major 
mergers of gas-rich galaxies, secular evolution (which includes both 
internal secular evolution and external secular evolution, the latter of 
which is driven by galaxy interactions), and hot halo accretion, which 
is presumed to be the dominate BH growth mode for low-excitation 
radio-loud AGNs. 



ski|2004 1 are predicted to drive the gas down to x 10- 
100 pc. 

The physics of gas inflow can be explored in some de- 
tail using multi-scale smoothed-particle hydrodynami- 
cal simulations. For example, adopting this approach 
Hopkins & Quataert ( 2010) 1 traced the dynamics of gas 
inflow over a broad range of scales (from x 10 kpc 
down to < 1 pc), taking into account the effects of 
stars, star formation, stellar feedback, and the gravita- 
tional influence of the BH. Their simulations revealed a 
much broader range of gas morphologies than typically 
expected, including spirals, rings, clumps, bars within 
bars, and nuclear spirals. All of these processes were 
effective at driving the gas down to ^ 10-100 pc but 
the gas inflow was found to be neither smooth nor con- 
tinuous. Any individual process was also found to be 
comparatively short lived, complicating attempts of re- 
lating morphological features to the presence of AGN 
activity, particularly on > 100 pc scales. 

Indeed, there is at best a marginal relationship be- 
tween the presence of AGN activity and galaxy bars in 
the local Universe, with AGNs and galaxies not host- 
ing AGN activity equally Ukely to host kpc-scale bars 



(e.g., 'Mulchaey & Regan 


1997; 'Knapen et al.| 2000 


|Martini et al. 


2003 


Simoes Lopes et al. , 2007]l; how- 



ever, we note that a relationship is found for narrow-Une 
Seyfert 1 galaxies (see §3.2.3). On smaller ^ 100 pc 
scales, high-resolution imaging of nearby systems has 
also failed to reveal any strong differences between 
the circumnuclear morphologies of AGNs and late-type 
galaxies not hosting AGN activity, with nuclear bars, 
spirals, discs, and circumnuclear dust equally likely to 



,Malkan, ^2004]). 



Pogge & Martini 2002 Martini et al. 2003) H 



2002 



unt & 



By contrast, a difference in the cir- 
cumnuclear regions between AGNs hosted in early type 
galaxies and inactive early type galaxies is found, with 



the AGNs only found in dusty environments (Simoes 



Lopes et al. 2007| l, suggesting that dusty structures are 



a necessary but not sufficient condition for BH growth. 
Although imaging studies have not revealed a morpho- 
logical signature on 100-pc scales that uniquely reveals 
the trigger of AGN activity, triggers of AGN activity and 
the identification of significant gas inflow may be more 
discernable using spectroscopic data. 

Integral Field Unit (IFU) observations provide spa- 
tially resolved spectroscopy and are a key tool to deter- 
mine the dynamics of gas and stars in nearby systems. 
IFU observations of nearby AGNs have revealed clear 
differences between the gas and stellar kinematics on 
< kpc scales, with the gas velocity fields showing signif- 
icant departures from the regular circular orbits traced 



by th e stars (e.g.,|Dumas et al.)|2007 Stoklasova et al. 
2009 Rodrfguez-Zaurfn et al. 2011| l. The gas kinemat- 



ics uncover a broad range of features not revealed by 
the photometric data, including warps, counter-rotating 
discs, and motions out of the galactic plane, likely due 
to gas inflow and outflow. However, these features are 
also seen in galaxies not hosting AGN activity, indicat- 
ing that the gas inflow that feeds the BH must occur on 
yet smaller scales. 

IFU observations on a: 10 pc scales of several nearby 
AGNs are starting to reveal nuclear gas inflows, when 
contributions from the stellar potential are carefully 
modelled and removed from the ionised gas kinemat- 



ics (e.g., Fathi et al. 2006[|Dumas et al. 2007 Storchi- 


Bergmann et al. 


'2007 i'Riff^el etal. 


'2008 Daviesetal. 


,2009, SchnoiT 


Miiller et al. 


2011). The strongest 


evidence for a 


nuclear gas 


inflow to date is from 



NGC 1097, where the gas is seen streaming along spi- 
ral structures towards the nucleus on ^ 3.5 pc s cales, 

60 km s"^ (e.g., Fathi 



with inflow velocities of up to 

'et al.| |2006[ [Davies et aL][2009l see also [Prieto et al.' 



|2005i. However, a significant restriction in the identi 
fication of nuclear-scale gas inflows is that the gas kine 
matics of luminous AGNs are often dominated by gas 
outflow signatures (e.g., Storchi-Bergmann et al.|20T0 



Hicks et al. in prep), and therefore the best constraints to 
date are limited to low-luminosity AGNs where the gas- 
inflow rates are comparatively low. This observational 
challenge can be overcome, in at least some cases, us- 
ing spectropolarimetry, which can help distinguish be- 
tween structures with different geometries (such as gas 
inflow from gas outflow) even in spatially unresolved 
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regions (e.g., Young 2000 Smith et al. 2004[ l. For ex- 
ample, using spectropolarimetry Young et al. (in prep) 
identified small-scale gas inflows in the Seyfert 1 galax- 
ies NGC 4151 and Mrk 509, with inward radial veloc- 
ities of a; 900 km s"' and mass inflow rates of x 0.2- 
0.7Moyr-i. 

IFU observations have revealed the signatures of gas 
inflow down to ^ 10 pc scales in several nearby AGNs. 
However, paradoxically, the gas inflow rates estimated 
from the ionised gas are often significantly lower than 
the BH accretion rate (e.g., (Fathi et al. 2006 ; Storchi 
Bergmann etiTIIIOOTl [Riffel et al.||2008, , Davies e~ 



velocity dispersion (osph) as 
GMbh 



r < 



(2) 



>rchi-| 
5taiJ 



The density of the gas on these scales can be high 
and if it becomes unstable then it will collapse and form 
stars (e.g., [Toomre 1964 Kennicutt 1989 Thompson 
et al.j |2005| l, potentially robbing the BH of the fuel fliat 
it needs to grow. For example, an isothermal gas disc is 
expected to become unstable when 



2009 Schnorr Miiller et al. 2011 1. These discrepancies 
may indicate that the ionised gas provides a poor tracer 
of the inflowing gas rate, which is expected to be dom- 



(3) 



which is equivalent to Q < 1 in terms of the Toomre 



inated by cold molecular gas (e.g. Davies et al. |2009 



[Hicks et al.| [20091 IMiiller Sanchez et al.| |2009| l; the 
ionised gas may be the outer skin of the molecular gas 
inflow (e.g., Riffel et al. 2008[ l. Indeed, the mass inflow 
rates estimated using molecular gas tracers in several 
nearby AGNs are large enough to power the observed 



parameter (e.g., Toomre[ 


1964[ [Thompson et al.|[2005| 


McKee & Ostriker 2007 


1. In this equation Cj is the ef- 



fective sound speed in the gas, k is the epicyclic fre- 
quency, and 2g is the gas surface density. 

These size scales correspond to that expected for (and 



in some cases observed; e.g., Tristram et al.|2007 2009 



et al. 2009| l. The molecular gas luminosities of AGNs 



AGN activity (e.g., [Davies et alT| [2009{ [Miiller Sanchez] [Raban et aL][2009[ [Kishimoto et al.[[2011| l, the dusty 



on < 100 pc scales are also potentially larger than those 
found in galaxies not hosting AGN activity, suggesting 
a greater reservoir of cold circumnuclear gas in at least 



some AGNs (e.g., [Miiller Sanchez et aL][2009[ Hicks 
et al., in prep). 

The current concensus view is therefore one where 
gas inflow from kpc scales down to the central ^ 10- 
100 pc region occurs in all gas-rich galaxies. The gas is 
driven into the central region through a series of gravita- 
tional instabilities, which helps to overcome the angular 
momentum of the gas. The clearest triggers of AGN ac- 
tivity in the local Universe are on the smallest scales, 
which is perhaps not surprising given that the gas in- 
flow time from the host galaxy down to the pc-scale of 
the BH accretion disc is likely to be i 10^ years, com- 
parable to or longer than the duration of the AGN phase 



(e.g.. Martini et al. 2003 Hopkins & Quataert 2010 1 



2.2. The competition between star formation and black- 
hole accretion 

On scales of order » 1-10 pc the inflowing gas is ex- 
pected to hit a bottleneck. The axisymmetric structures 
that have driven the gas from ^ 10 kpc down to ^ 10 pc 
are likely to have a negligible effect on the gas in these 



molecular torus that is thought to obscure the central 
regions of the AGN in inclined systems (e.g., as postu- 
lated in the unified AGN model; [An"tonucci| 1 993) . That 
these structures should be on similar scales is unlikely 
to be a coincidence - the dense gaseous and dusty struc- 
ture is potentially both the outer regions of the BH fuel 
supply and a stellar nursery. Indeed, the change from the 
galaxy potential to the BH potential can produce warps 
in the gas that could obscure the accretion disc (e.g., 
Hopkins & Quataert' '20T0'i 'Lawrenc e & Elvisj [2nT0| ); 
see also Bregman & Alexander (2009 1 for warping of 
the accretion disc due to dense cusps around the BH. 
However, even at these small scales, the gas still has to 
lose ^ 99% of its angular momentum before it reaches 
the sub-pc scale of the accretion disc. 

Numerous studies have shown that young stars are of- 
ten found in the central regions of nearby AGNs (e.g., 
Storchi-Bergmann et al.[ [200T] [2005] [Cid Fernandesl 
et al. 2004| l, suggesting a connection (causal or non 
causal) between AGN activity and star formation. Using 
spatially resolved spectroscopy of the central regions 
(ss 10-100 pc) of nine nearby Seyfert galaxies, Davies 



et al. ( 2007) provided compelling evidence for a causal 
cormection between AGN and star-formation activity in 



inner regions (e.g., Hopkins & Quataert 2010 1. It is also 



on these scales that the gravitational influence from the 
BH will overcome that of the host galaxy; this radius is 
determined from the BH mass (Mbh) and the spheroid 



at least some systems (see also Gonzalez Delgado et al. 
2001) ■ [Davies et aL] ( [2007| l showed that flie peak AGN 
activity occurs =s 50-200 Myrs after the onset of star 
formation (see also 'Schawi nski et al.|2009[ [Wild et aL] 
2010, and §3.2.3). The time delay between the star for- 
mation and AGN activity suggests that the star-forming 
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region may provide the fuel for the BH, via winds from 
massive stars and supernovae (e.g., VoUmer et al.|2008 
however, see |Hopkins|20 1 1 | l . 

These studies imply a complex interplay between 
AGN activity, star formation, and stellar winds in the 
vicinity of the BH. A wide range of analytical mod- 
els, numerical simulations, and hydrodynamical simu- 
lations indeed predict that stellar winds and supernovae 
will enhance the BH mass accretion by injecting turbu- 



lence into the gas disc (e.g., |Wada & Norman 2002 



ISchartmann eTaTl [2009) |Hobbs et al.| |2011| l. How 



ever, the magnitude of the effect depends on the as- 
sumed prescriptions of the stellar-related outflows and 
the gas accretion timescales. For example, if the gas- 
accretion timescale is long then the majority of the gas 
will collapse and form stars rather than be directly ac- 



creted onto the BH (e.g.. King & Pringle 2007 Nayak- 



shin et al. 2009 | l. In this scenario the gas accreted onto 



the BH will be predominantly from stellar-related out- 
flows (e.g., SNe; AGB stars) and is therefore triggered 
after the initial starburst activity. Assuming that the 
stellar-mass loss is dominated by AGB stars, the BH 
accretion-star formation ratio would be « 10^^-10"^ 
(e.g., |Jungwiert et aL] 12001) 1. By contrast, on the basis of 



their multi-scale hydrodynamical simulation, Hopkins 



& Quataert ( 2010) l argue that a precessing eccentric disc 
will form within ^ 10 pc of the BH in gas-rich systems, 
which can efficiently deliver gas to the accretion disc. 
The rate of BH accretion and star formation within the 
central 10 pc in the [Hopkins & Quataert| (2010i model 
is approximately equal, and therefore 2-3 orders of 
magnitude larger than that predicted from stellar-mass 
loss. 

While there is disagreement over how the gas is deliv- 
ered to the accretion disc from « 10 pc, it is clear that the 
gas inflow is complex and cannot be easily prescribed 



as a simple equation (e.g., Hopkins & Quataert 2010 
Power et al. 201 1). A popular approach in the treat- 
ment of BH accretion in galaxy formation simulations is 
the Bondi-Hoyle mefliod ( [Bondi & Ho"yie)[T944)[Bondi| 



1952) 1, which assumes spherical accretion onto the BH 



and therefore doesn't consider angular- momentum lim- 



itations (e^g^jDT^atteo^Firi [2005 1 [Springel et al. 



2005a Boofli & Schaye 'Mjg). To address this sim 



plification, Power et al. (2011 ) have produced a phys- 
ically self-consistent BH accretion model suitable for 
galaxy formation simulations that accounts for the an- 
gular momentum of the gas. The difference between 
their model and Bondi-Hoyle accretion is quite striking 
- with the Bondi-Hoyle method, the BH is accreting al- 
most continuously in gas-rich systems at high mass ac- 
cretion rates while, by contrast, a sizeable fraction of 



the gas in the Power et al. (2011 1 model collapses to 
form stars and the mass accretion rate onto the BH is 
comparatively negligible. 

The effect of stellar-mass loss on the BH accretion 
can be explored in even greater detail using small-scale 
hydrodynamical simulations of a nuclear starburst in 
the vicinity of the BH. The current suite of simulations 
(e.g., [Wada elaL) [2009) [Schartmann et aL) [2009) pOTO) l 
generically predict that stellar and supernova ejecta in- 
ject turbulence into the accretion disc, "puffing" it up 
and forming an optically and geometrically thick molec- 
ular torus - a single supernova can effect the whole torus 
and enhance the mass-accretion rate ( [Wada et ar)[2009| l. 
Radiative heating and pressure can slow down or stop 
the accretion process while radiative drag can enhance 
accretion. [Kawakatu & Wada[ ( |2008) l show that the frac- 
tional amount accreted decreases with an increase in the 
fuel supply since star formation in the disc itself can 
dominate over the growth of the BH; indeed, the pres- 
ence of nuclear stellar rings close to the BH in nearby 
quiescent galaxies provide indirect evidence in support 



of this hypothesis (e.g., Nayakshin & Sunyaev 2005 



Bender et al. , 2005 ; Paumard et al.[ |2006[|Schodel et al. 
,2007, Hopkins & Quataert, 201 0| l. 

The hydrodynamical simulations also predict that the 
gas can cause significant time-variable obscuration in 
the vicinity of the BH (e.g., [Wada & Norman, ,2002) 
[Wada et al.|[2009) see [Hicks et al.||2009| for observa- 
tional constraints). The predicted column densities and 
time-variable obscuration are in good agreement with 
those measured in nearby obscured AGNs (A^h > 10^^- 
10^^ cm"^ with absorbing column changes on < 1 year 



timescales; e.g., ,Risaliti et al. 1999 2002 Guainazzi 



et al.[|2005 [Cappi et al.|[2006 1. Taking the properties 



derived for the Seyfert 2 galaxy NGC 1068 from [Davies[ 
etaL[ ( [2007] l, [Schartmann et al.[ ( [20T0| used a hydrody- 



namical simulation to make detailed predictions of the 



properties of the obscuring torus (see Table 1 of Schart- 



mann et al.|2010 for the initial parameters). They found 



that the torus can be built up through the ejection of en- 
ergy from stellar winds. Their model reproduced the 
dense obscuring structure observed from interferomet- 
ric mid-IR observations on a: 0.5-1 pc and compares 
well with the extent and mass of the rotating disc in- 



ferred from H2O maser observations (e.g., GreenhiU & 
|Gwinn|[T997)|Raban et al.|[2009) l. 

Connections between the accretion onto the BH and 
obscuration in AGNs can be explored using X-ray data 
and mid-IR observations, which will probe the cen- 
tral source and the obscuring region, respectively. Us- 
ing high-spatial mid-IR imaging (typically probing < 



100 pc scales), Horst et al. (2008 2009 1 and Gandhi 
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et al. ( 2009 1 found a tight correlation between the unre- 
solved mid-IR core and the absorption-corrected X-ray 
luminosities of nearby AGNs, even in the presence of 



extreme Compton-thick (A^h > 10 cm ) absorption 
(see |Lutz et al.|200 Tfor similar constraints from mid-IR 
spectroscopy). These results indicate that the obscur- 
ing structure is directly heated by the AGN and further- 
more they suggest that the covering factor of the ob- 
scuring material is broadly constant (at a given AGN lu- 
minosity; see Lawrenceje.g., 1991 Simpson e.g., 2005 



for evidence of a luminosity-dependent covering fac- 
tor). These results are inconsistent with the simplest 
dusty torus models, which assume a smooth distribution 
of obscuring dust (e.g.. Pier & Krolik 1992; GranatoJ 
|& Danese|[T994l|Efstathiou & Rowan-Robinson,, 1995j ) 
and therefore predict that the mid-IR emission will be 
strongly absorbed in systems when the torus is inclined 
to the line of sight. However, these results are con- 



sistent with clumpy torus models (e.g., Nenkova et al. 



2002; Elitzur & Shlosman, 2006; Nenkova et al., 2008 
Schartmann et al. 2008; Honig et al., 2006; Honig &■ 

^^^^^^_l_l^_^_J ■ i I S i ■ I I ^zs I 

Kishimoto 2010[l, which predict no inclination-angle 



dependence on the strength of the mid-IR emission (see 
[Mullaney et"aL]|20 1 1 a| for full IR SED constraints that 
also support this view); dust in the narrow-line region 



can also lead to isotropic dust emission (e.g., Efstathiou 
[etaLl [1995, Schweitzer et al.|[2008) . 

The results presented in this section show that it is 
challenging to accurately predict the fate of the gas on 
=s 1-10 pc scales. Generically, the current models pre- 
dict that the gas inflow on these scales will lead to a 
complex interplay between AGN activity, star forma- 
tion, and stellar winds. However, different models pre- 
dict different rates of mass accretion and star formation. 
The current view of the AGN dusty torus is that of a 
structure more dynamic than that originally postulated 
in the unified AGN model (e.g., Antonucci 1993| l. The 
AGN "torus" is potentially both the outer gas reservoir 
for the BH accretion disc and a nuclear stellar nursery, 
built from the gas inflow and stellar-related winds and 
outflows. 

2.3. Mass accretion onto the black hole 

The power house behind AGN activity is mass accre- 
tion onto the BH on <K 1 pc scales. The basic theory of 
the accretion disc (the so-called a disc) has been around 



for almost 40 years (e.g., Shakura & Sunyaev 1973 1. In 



the case of optically thick accretion, viscosity (the a pa- 
rameter, which is a major cause of uncertainty) causes 
the gas to lose angular momentum and fall towards the 
BH, transporting the angular momentum outwards. The 
optically thick accretion disc is geometrically thin and 



the energy spectrum is comprised of multi-temperature 
components, which reach their peak temperature at the 
centre (e.g., |Pringle][T98T][Rees|ri984) . The local efi'ec- 
tive temperature approximates to (e.g., Goodman|2003[ l 



Teff^ 6.2x10-1 — - K. 



(4) 



where A is the Eddington ratio (Lboi/^^Edd; defined be- 
low), eo.i is the radiative efficiency in units of 0.1, Mg 
is Mbh in units of 10*^ M©, r is the radius, and 7?s is the 
Schwarzschild radius (i.e., the event horizon), defined 
for a spherical non-rotating BH as 



2GM, 



BH 



(5) 



The Eddington ratio (A - Lboi/i'Edd) indicates the rela- 
tive growth rate of the BH based on the AGN bolometric 
luminosity (Lboi) and the Eddington luminosity (Z^Edd), 
defined as 



Z'Edd - 



AnGMBunipC 



(6) 



where nij, is the mass of the proton and crj is the cross 
section of the electron. The Eddington luminosity is 
achieved when the outward radiation pressure equals the 
inwards gravitational force - if the Eddington luminos- 
ity is exceeded then (in the absence of an additional in- 
ward pressure) the gas is expelled; see Fig. 1 for the 
Eddington ratios of some nearby AGNs. 

An alternative solution to the accretion-disc equa- 
tions is the optically thin accretion disc, a state achieved 
at low mass accretion rates and generically referred to 
as Radiative ly Inefficient Accretion Flows (RIAFs; e.g, 
[Narayan & Yi|1994[pandford & Begelman'1999). The 
optically thin accretion disc is unable to cool efficiently 
(the cooUng time exceeds the accretion time) and the 
energy is lost through non-radiative processes such as 
convective transport of energy and angular momentum 
to large radii or an outflowing wind. Magnetic fields are 
likely to dominate viscosity and the transport of angular 
momentum in accretion discs, as suggested by numeri- 
cal simulations of Magneto-Rotational InstabiUty (MRI; 



e.g., |Balbus & Hawley|199T][T998HBaibus|2003| ). 

Depending on the accretion disc viscosity and the 
mass accretion rate, the accretion disc is predicted to 
change between an optically thick and an optically thin 
state - the transition from optically thin to optically 
thick is likely to co rrespond to an Eddington ratio of 
-10-2 (e.g. 



10 



-3 



Esin et al. 



1997 



Gallo et al. 



2003 



[Maccarone) |2003| l. However, despite the accretion disc 
being the ultimate source of power in AGNs, we lack the 
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observational data to robustly test these models. There 
are hints that the emission from AGNs deviates from the 
simplest accretion-disc models. For example, micro- 
lensing of the central parsec region of the accretion 
disc in a number of quasars (e.g., Pooley et al. 2007 



Floyd et al. 2009^ have shown that the continuum emis- 
sion profile drops more quickly that than expected from 
the basic a accretion-disc model but is broadly consis- 
tent with that predicted from MRl accretion-disc models 



(e.g., |Agol & Krolik|[2000 i. 

Fortunately, significantly better tests of the accretion- 
disc models are available for accreting stellar-mass 
BHs, so called X-ray binaries (XRBs), for several key 
reasons: (1) the data is typically of a significantly higher 
signal-to-noise ratio than that found for AGNs since 
many XRBs lie in the Galaxy, (2) the accretion-disc 
state of XRBs can change on timescales of less than a 
day in response to changes in the mass accretion rate, 
allowing for different accretion-disc components to be 
identified and analysed; by comparison the accretion- 
disc state of AGNs will vary on timescales of centuries 
(the timescale is a function of the BH mass; [Mirabel &] 
Rodriguez ][l998')0and (3) the accretion disc of XRBs 
peaks in the X-ray band (» 1 keV; a: 10^ K; see Eqn. 4), 
where it can be studied in great detail using the cur- 
rent generation of X-ray telescopes. By comparison the 
accretion disc of AGNs peaks at far-UV wavelengths 
(a 10 eV; X! 10^ K; see Eqn. 4), where the emission 
is absorbed by the interstellar medium in the Galaxy. 
Despite these observational challenges in comparing 
AGNs and XRBs, there is good evidence that the accre- 
tion process is essentially the same for both small and 
large BHs and we can therefore use studies of XRBs to 
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Figure 3: Schematic diagrams to illustrate the expected range of BH 
accretion modes as a function of Eddington ratio. The orange-red re- 
gions refer to the optically thick accretion disc and the blue regions 
refers to the X-ray "corona". Adapted from the figure produced for 
X-ray binaries by Done, Gierlinski, & Kubota (2007) and revised to 
show how the various optical spectroscopic classes of AGNs might 
correspond to the different BH accretion modes and the dominate re- 
lease of energy (radiative or mechanical; i.e., predominantly due to 
jets). 



accretion rate (e.g., Nowak 1995 Done et al. 2007| l. 



better understand mass accretion in AGNs (e.g., perH [Done et al .1120071 for a more in-depth discussion). 



loni et al.| [2003| [Falcke eFaLl [20041 [McHardy et al 
2006| 



At high Eddington ratios (A i 10^) the X-ray emission 
is often dominated by a disc-like spectrum, assumed 
to be due to an optically thick accretion disc, with a 
weaker power-law continuum component. At low Ed- 
dington ratios (A t 10~^) the hard X-ray spectrum is 
dominated by the power-law continuum and has a weak 
disc-like spectrum; however, the conversion from the 
high-soft to the low-hard state is more complex than 
that briefly implied here and should not be taken as 
definitive (e.g., see Remillard & McClintock[[2006 and 

In 



2.3.1. Accretion onto stellar-mass black holes 

The X-ray spectra of XRBs are typically classified 
into "X-ray states", which qualitatively describe their 
behaviour and broadly relate to the properties of the ac- 



cretion disc (e.g., see Table 2 of Remillard & McClin- 



tock|2006l l - a low-hard state, where the majority of the 
X-ray emission is produced at > 10 keV and a high- 
soft state, where the emission is strong at < 10 keV but 
weak at > 10 keV. An individual XRB can be seen to 
traverse from one X-ray state to another on timescales 
of less than a day in response to changes to the BH 



the high-soft state the accretion disc is optically thick, 
and therefore viscous, and emits a thermal spectrum, 
while in the low-hard state the accretion disc is opti- 
cally thin and radiatively inefficient and the X-ray emis- 
sion is dominated by the inverse Compton scattering of 
photons from the outer regions of the accretion disc (but 
can also have significant contributions from synchrotron 
emission and synchrotron self-Compton emission). Fur- 
ther X-ray states have also been identified, such an in- 



termediate very-high state (see Table 2 of 'Remillard & 



[McClintock|2006 ]), which likely correspond to differing 
contributions from the accretion disc, inverse Compton 
scattering, and an associated jet. 

The radio emission from XRBs is found to be linked 
to the X-ray state (e.g.. Fender et al. 20041 1. The radio 



'The X-ray variability of AGN seen on intra-day timescales is re- 
lated to the mass accretion but is unrelated to the accretion-disc state 
(e.g., [Uhich"erari[T997l[Vaughan et al"][2003l[Uttley et al.[[2005) . 



emission is bright in the low-hard state and increases 
in luminosity in response to an increase in the Edding- 
ton ratio up until the low-hard-high-soft state transi- 



tion at /i a: 10 ^. The radio emission dramatically 
declines when the XRB traverses to the high-soft state 
(e.g., I Gallo et"aL] [2003] [Fender et al.| [2004| . The radio 
emission is consistent with that expected from a syn- 
chrotron emitting jet and the change in the radio power 
is possibly due to the scale height of the accretion disc 
(which is geometrically thick in the optically thin state) 
and provides the seed electrons to drive the jetj^ The jets 
in XRBs generate large amounts of mechanical energy 
which can have a significant effect on the surrounding 
interstellar medium. 

The most luminous XRBs have X-ray luminosities > 
10^' erg s ' and are called ultra-luminous X-ray sources 
(ULXs). ULXs hkely comprise a heterogenous popu- 
lation (se e [King et al.||2001[ |Begelman|[2002| [Roberts 



[2007] and [Gladstone et al.||2009l for a broader discus 
sion of the nature of ULXs), of which a sizeable fraction 
are stellar-mass BHs accreting at super-Eddington rates. 
ULXs therefore provide direct insight into the proper- 
ties of accretion discs at extreme mass accretion rates. 
The X-ray spectra of many ULXs can be fitted with an 
accretion-disc spectrum and an optically thick "corona" 
(e.g., Stobbart et al. 2006). However, the identification 
of the soft X-ray emission from the accretion disc in the 
presence of the large amounts of absorption inferred by 
the optically thick corona appears to be contradictory. 
A more physically plausible explanation is that the op- 
tically thick corona is an outflowing wind, launched in 
the vicinity of the accretion disc but more extended than 
the accretion disc and therefore not directly obscured. 
Hydrodynamical simulations of BHs accreting at super- 
Eddington rates indeed predict the presence of a radia- 



tively driven high-column density wind (Ohsuga et al. 



2005 Ohsuga 2007 1, providing indirect support for this 



model. The potential impact of these winds can be seen 
from the production of "bubbles" or emission-line "neb- 
ulae" over * 100 pc scales (e.g., Roberts et al. 2003 
ikuU et al. 2010| l. These emission-line nebulae may 
be analogous to those seen in AGNs (i.e., the narrow- 
line region; e.g., [Wilson et al.|1988| [Schmitt & Kinney| 
[T9951 l. 

2.3.2. Accretion onto massive black holes 

The accretion-disc in AGNs peak at far-UV wave- 
lengths, where we lack direct observational constraints 
due to absorption by interstellar gas in the Galaxy. How- 
ever, the properties of the AGN accretion disc can be 



Strong radio-emitting XRBs have traditionally been called micro- 
quasars (e.g., Mirabel et al. 1992i |Mirabel & R odriguez 1994 1998) 
but it now appears likely that they are just the radio-bright subset of 
the XRB population. 



indirectly inferred from the strength of the UV and op- 
tical emission lines. The UV continuum is efficient at 
photoionising the gas since the majority of the atomic 
transitions occur at energies corresponding to UV wave- 
lengths. On the basis of this, the low-ionisation emis- 
sion lines produced by Low-Ionised Nuclear Emission- 
line Regions (LlNERs; e.g., Heckman 1980 Baldwin 



[et al.[[T98 T) would correspond to the low-hard power- 
law dominated state while the high-ionisation emis- 
sion lines produced by Seyfert galaxies and quasars 
would correspond to the high-soft accretion-disc domi- 



nated state (e.g., Kewley et al. 2006 Sobolewska et al. 



[20TT] i; see Fig. 3. 

A major source of uncertainty in AGN studies is the 
mass-accretion rate, which requires an accurate estimate 
of the AGN bolometric luminosity. Since the majority 
of the accretion-disc emission in AGNs is produced at 
far-UV wavelengths, the AGN bolometric luminosity is 
typically determined using the continuum emission in 
a given wave band (L^; typically at X-ray or optical 
wavelengths) with a wavelength dependent bolometric 
correction factor (BCa) applied to calculate the total ra- 



diative output (e.g., Elvis et al. 
|2004,, Hopkins et aL|[2007] i; 



1994 Marconi et al. 



I.e., 



^bol - L^BCa- 



(7) 



The majority of studies assume an optically thick ac- 
cretion disc when determining BC. However, as shown 
for XRBs, an optically thick accretion disc is only 
likely to be valid at comparatively high Eddington ra- 
tios (A i 10^^-10"^) and therefore the mass-accretion 
rates estimated in studies of low-Eddington ratio AGNs 
will be wrong. Indeed by fitting the X-ray-optical SEDs 
of nearby AGNs, interpolating at UV wavelengths, ' Va-[ 
[sudevan & Fabian (2007 , 2009 ) have shown that the 
AGN bolometric correction may be a strong function 
of the Eddington ratio: the X-ray-bolometric conver- 
sion factors change from ^ 15-60 for Eddington ratios 
of A X 10^^-1. The dominant factor in variations in the 
SEDs and emission-line properties of quasars may also 



2000 



be due to the Eddington ratio (e.g., |Sulentic et al. 
[Kuraszkiewicz et al.|[2009} . 

Even greater insight into the bolometric output of 
AGNs can be made by carefully fitting the X-ray and 
optical spectra with physically motivated models. Re- 
cent work has focused on fitting the X-ray spectra of 
AGNs with an accretion disc, power law, and Compton 
component while taking account of the various stellar 
and AGN components in the optical spectra to search for 
connections between the underlying physical processes 



(e.g., Jin et al. 2009 Sobolewska et al. 2011 1. Hutton 
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et al. (in prep) have taken this approach for AGNs iden 
tified in the SDSS (e.g., [York et al.| |2000| [Abazajian 
|et al.j |2009| ) with excellent-quality X-ray data obtained 
from the XMM-Newton serendipity survey (e.g., Watson 



et al. 2009 1. A wide diversity of X-ray-optical SEDs 
are found, which would otherwise not be apparent from 
the X-ray data alone. There is some evidence for cor- 
relations of the strength of the X-ray emitting compo- 
nents and the width of the Wfi emission line. However, 
since AGNs vary significantly in the X-ray band, con- 
temporous X-ray and UV-optical data are required to ac- 



count for intrinsic variations (e.g., Vasudevan & Fabian 
|2009l l. 

2.3.3. Sgr A*: the closest accreting massive black hole 
The closest massive BH to us is Sgr A*, which lies 
at a distance of » 8.3 kpc in the centre of the Galaxy 



(e.g., Reid 1993 Melia & Falcke 2001 Gillessen etal. 



|2009 1. The BH mass of Sgr A* has been measured 
to unprecedented accuracy using ^ 16 yrs of monitor- 
ing of the orbits of stars in the vicinity of Sgr A*: 



see 



Mb h = (4.3 + 0.2) X 1 0^ Mo ( |Gillessen et al.|2009 
also |Ghe z et al. 2008). Sgr A* produces emission over 
a broad waveband (radio-X-ray; e.g., |Marrone et al.| 
[20081 |Dodds-Eden"eraL][20lT| However, with an es 
timated Eddington ratio of A ^ 10"^-10^° (e.g. 



Melia 



& Falcke 2001 Eckart et al. 2006 1, Sgr A* is ^ 5-8 



orders of magnitude below the optically thick-optically 
thin accretion state and is only identifiable as "active" 
due to its proximity. For example, none of the AGNs 
and LINERs investi gated in the nearby galaxy sample 
of ( |Ho et al.| '1997a) have Eddington ratios of ^ < 10"* 
(see Fig. 9 of Ho 2008) , indicating that Sgr A* is in 
a quiescent state. Sgr A* therefore provides unique in- 
sight into low mass-accretion events onto a massive BH. 
The broad-band emission from Sgr A* appears to 



have at least two components (e.g., Dodds-Eden et al. 
201 l| l - a continuously emitting "quiescent" state and 
a flaring state. The flares are often quasi-periodic 



and can vary on short timescales (e.g,. IBaganoff" et al. 
[200T1 IGenzel et al. , IMT. Ghez et al.| |2004[ [Marrone] 



[eFaLl [2008; Dodds-Eden et al. 20lT); the radio-near- 
IR spectral energy distribution of the flares observed in 
Sgr A* are consistent with synchrotron emission from 



transiently heated electrons (e.g., [Markoff et al. 2001 


Gillessen et al. 


2006 Trippe et al. [2007| [Eckart et al. 


2009 Dodds-Ed 


en etal. 2010 201 l)l. The origin of the 



flaring emission is a topic of hot debate and may be due 
to an orbiting hot spot, magnetic heating/reconnection, 
stochastic acceleration processes, or accretion-rate en- 



hancements (e.g., Markoff^ et al.| 


2001 [Liu & MeUal 


|2002||Yuan et al.||2003||Dodds-Ed 


en etal.l 20111. The 



modest BH accretion rates inferred by the emission 
from Sgr A* are consistent with that expected from the 
mass loss of stars in the Galactic centre (e.g., |Cuadra[ 
[etaLl[2006l[2008| ). 

On the basis of the current mass accretion rate onto 
Sgr A*, it would take « 5-6 orders of magnitude longer 
than the Hubble time to grow the BH to its current mass, 
indicating that Sgr A* must have been orders of mag- 
nitude more active in the past. The identification of re- 
flected X-ray emission off" molecular clouds in the vicin- 
ity of Sgr A* provides a potential "X-ray echo" of past 
activity from Sgr A*. These data suggest that Sgr A* 
may have been up-to ^ 3 orders of magnitude brighter 
100 years ago (e.g., jMuno et al. 2007 [ Ponti et al. 



2010 1. However, these implied mass accretion rates are 



still extremely modest and Sgr A* must have been much 



more luminous in the past (Nayakshin et al. 2007 1; see 
[Micic et aL] ( |201 l| l for potential BH growth histories for 
Sgr A*. 



3. What are the links between black-hole growth 
and their host galaxies and large-scale environ- 
ments? 

In the previous section we discussed the mechanics 
and processes of gas inflow from the host galaxy down 
to the central BH, across > 5 orders of magnitude in 
size scale. The BH-spheroid mass relationship implies 
that this, otherwise seemingly unlikely journey of gas 
inflow, has occured in all spheroid-hosting galaxies at 
some point over the last ^13 Gyrs. The tightness of the 
BH-spheroid mass relationship farther suggests a con- 
nection between AGN activity and star formation over 
kpc scales. In this section we explore the triggering 
mechanisms and sites of AGN activity and investigate 
the connection between AGN activity, star formation, 
and large-scale environment. We also investigate the 
growth of BHs across cosmic time and explore how the 
properties and triggering of distant AGN activity differs 
from that found in the local Universe. 

3.1. The identification of AGN activity 

To accurately explore when, where, and how the 
growth of BHs has been triggered requires the identifi- 
cation of AGN activity across all environments and red- 
shifts. The presence of dust and gas in the vicinity of 
the BH and star-forming regions along the line of sight 
means that penetrating observations are required to un- 
ambiguously reveal the signatures of AGN activity in all 
systems. The emission from star formation and starlight 
in the host galaxy can also dilute and mask the emission 
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from the AGN. These comphcations make the construc- 
tion of a complete census of AGN activity a significant 
challenge. 

The most extensive studies of AGN activity in the 
local Universe have utilised optical spectroscopy (e.g., 



can select AGNs almost irrespective of the presence of 



Huchra & Burg 1992 |Maiolino & Rieke 


1995 |Ho 


et all |1997a|c|b| Kauffmann et all |2003 


Heckman 



etaL] [20041 [Hoi |2008'). Since the AGN narrow-line re- 
gion extends beyond the obscured central region 10- 
1000 pc versus » 1-10 pc), optical spectroscopy is able 
to identify even heavily obscured AGNs, so long as the 
AGN narrow emission lines are not obscured by dust in 
the host galaxy. Sensitive radio, mid-lR spectroscopy, 
and X-ray observations (particularly at > 10 keV; e.g., 
jSazonov et an[2007l [Winter et al.|[2009l [BuHon et al.| 
[201 l| l further extend our understanding and census of 
nearby AGN activity, providing the identification of 
AGNs even in the presence of significant host-galaxy 



dust and star formation (e 




Ho et al. 


2001 Ulvestad 


& Ho, 2002; Filho et al. 


2006; Satyapal et al. 


2008'; 


Goulding & Alexander, ,2009 


,Zhang et al. 


,2009,, Win-, 


ter et al.| 2010| [Nardini et al.| 2010 1. 


From a com- 



bination of these multi-wavelength identification tech- 
niques, we are now close to a complete census of AGN 
activity in the local Universe down to a given sensitivity 
limit (i.e., lower-luminosity AGNs can still be unidenti- 
fied). However, even intrinsically luminous AGNs can 
remain unidentified if the obscuration towards the BH 
is extremely high and the central region is heavily ob- 
scured by dust in the host galaxy (see the conflicting ev- 
idence forA^^^220jClements et al.|2002[ Spoon et al. 



[20041 [Iwasawa et al.]2005t|Downes & Eckart 2007|). It 
is important to bear in mind these potential limitations. 

Optical spectroscopy is also effective at identifying 
distant obscured AGNs (e.g., [S teidel et al. , 2002'; 'Pol-' 
[letta et al.| [2006t [Alexander etal.,, 2008b,, Juneau et al.. 



2011 Yanetal. 2011 1. However, since many of the key 
AGN emission-line diagnostics move into the near-IR 
band at z i 0.4, optical spectroscopy alone can be quite 
limited in identifying large numbers of AGNs (but see 
lJuneau et al.|20TTl|Yan et al.|201 ll[TYouille et al.|201 1| 
for techniques that extend the utility of optical spec- 
troscopy for AGN identification out to z ~ l)j^ Cur- 
rently the most efficient and effective identification of 
distant AGNs is made with X-ray observations, which 



'The optical signatures of distant AGNs are also more easily di- 
luted from host-galaxy emission than nearby AGNs due to the larger 
angular-size distance for distant systems and the overall increase in 
star-formation activity at higher redshift ( e.g., [Noeske et aT [2007 



Daddi et al.||2007b[[Pannella et al.| [20091 [Rodighiero et al.|(20T0a 
Elbaz etal.||2011] 



obscuration (e.g., Bauer et al. 2004 Tozzi et al. 2006 



[Gilli et aL] [2007| i. X-ray observations are effective at 
identifying AGN activity because the X-ray emission 
from star formation is typically weak; the positive K 
correction at X-ray energies for distant obscured AGNs 
also means that < 10 keV observations are particularly 
effective at identifying heavily obscured AGNs at high 
redshift (e.g., [Tozzi et al. 2006 Alexander et al. 2011 



Comastri et al. 2011 Feruglio et al. 2011 1. Sensitive 
IR and radio observations can farther extend the cen- 
sus of distant AGN activity by identifying AGNs where 
the absorbing column densities are so high that not even 



X-ray photons can escape (e.g.. Lacy et al. 2004; Stern 


et al. 2005[ Donley et al.l 


2005, ,2008, ,Alonso-HeiTero 


et al.[ 20061 Daddi et al.[ 


2007a| Hickox et al.| 2007| 



,20lT] l. However, since star formation can also produce 
luminous IR and radio emission, it is often challeng- 
ing to disentangle the AGN emission from that of the 



host galaxy (e.g., PoUetta et al. 2007 MuUaney et al. 
|2011a^ . 

Differences in the various approaches in AGN iden- 
tification can also make it challenging to reconcile re- 
sults from different studies. It is therefore important to 
always consider how obscuration or host-galaxy emis- 
sion can effect the completeness of any AGN selection 
technique. In the following sub sections we explore the 
processes of AGN activity using observations at X-ray, 
optical, IR, and radio wavelengths. We investigate the 
ubiquity of AGN activity, the AGN triggering mecha- 
nisms, the connection between AGN activity and star 
formation, the evolution of AGN activity with redshift, 
and the role of environment in the triggering of AGN 
activity. 

3.2. AGN activity in the local Universe 

3.2.1. Where are the massive black holes growing? 

AGN activity is common in the local Universe, with 
~ 5-10% of nearby galaxies found to clearly host op- 
tical AGN activity (i.e., identified as Seyfert galax- 
ies from their optical emission-line properties; e.g., 
'Veifleux & Osterbrock'1987VKewley et al. '2001 '20061 
Maiolino & Rieke 1995, Ho et al. 1997a; Hao et'dl 



2005 The optical AGN fraction can be significantly 



higher if galaxies with Low-Ionisation Nuclear Emis- 



sion Regions (LINERs; [Heckman ^ 1980 ) are also taken 



into account (e.g.. Ho et al. 1997a). However, LIN 



ERs comprise a heterogenous mix of composite AGN- 
star-forming galaxies, optically obscured AGNs, low- 
luminosity AGNs, and early type galaxies without AGN 
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Figure 4: AGN fraction as a function of BH mass for different AGN 
selection approaclies: optically identified luminous AGNs (filled 
squares; Best et al. 2005), lower-luminosity AGNs (crosses; Gould- 
ing et al. 2010), and radio-luminous AGNs (open circles; Best et al. 
2005). 



activity (e.g.JBinette et al.||1994||Ho et al.||2001 



2003 



Eracleous et al.' '2002', 'Satyapal et all, 120041 IK ewley| 
et al. 2006 ; Shields et al. 2007 ; Gonzalez -Martin et al.l 



2009 Goulding & Alexander 2009 Sarzi et al. 2010 



Capetti & Baldi 2011 1, leaving some ambiguity over 
the true AGN fraction in the LINER population. Mid- 
IR spectroscopy can provide a more complete census 
of nearby AGN activity than that obtained from optical 
spectroscopy due to the identification of high-excitation 
emission lines (e.g., [Ne V]/114.3 fim) in systems where 
the narrow-line region is obscured by dust in the host 



galaxy (e.g., | Satyapal et al. 2008 Gouldmg & Alexan 



exan-| 



|der|[2009HTommasin et al.r 20lO; Petric et al. 201 
However, at present, the number statistics obtained from 
high-resolution mid-lR spectroscopy are poor compared 
to those achieved from optical spectroscopy (samples of 
tens rather than thousands of objects). 

The fraction of galaxies hosting AGN activity pro- 
vides a measure of the duty cycle of BH growth. The 
fraction of local galaxies hosting AGN activity as a 
function of BH mass is shown in Fig. 4. The frac- 
tion of galaxies hosting optically luminous AGN activ- 
ity (L[oiii] > 3 X 10^ Lq, equivalent to a 2-10 keV lu- 
minosity of Lx i lO'*-' erg s"', assuming the conversion 
of Mulchaey et al.|1994[ i is approximately flat at » 5% 
over Mbh ~ (0.3-10) xlO^ Mq and drops substantially 
to lower and higher BH masses (e.g., |Best et al.j [2005[ ). 
The AGN fraction rises towards lower AGN luminosi- 



ties, as expected since there will always be more low- 
accretion rate systems than high-accretion rate systems. 
For example, the fraction of galaxies hosting AGN ac- 
tivity =s 3 orders of magnitude fainter than the opti- 
cally luminous AGNs is =a 20-30% over Mbh ~ (0.1- 



10) xlO^ Mq (e.g., [Goulding et al.| pMO} . Note that 
Sgr A*, the active BH at the centre of the Galaxy is still 
Si 4 orders of magnitude fainter than these sensitivity 
limits! By contrast, the fraction of galaxies hosting lu- 
minous radio AGN activity (^i.4ghz 

> 1024 W Hz-1) 

is strongly dependent on the BH mass, rising from 
i 0.1 % at Mbh < 10*^ Mq (where optical and mid- 
IR identified AGNs are comm on) to « 10% at Mbh > 
lO'' Mq (e.g., iBestetaT, "2005). Luminous radio AGNs 
represent the minority of the AGN population at all red- 
shifts (radio-loud AGNs account for t 10% of the total 
AGN energy budget and cosmological BH growth) but 
they are likely to have played a crucial role in the evolu- 
tion of the most massive galaxies and BHs (e.g.,|Bower| 



let al.|2006 l |Croton et al.|2006l [Merloni & Heinz|2008 



see §5.1)|^ 

These suite of results show that mass accretion onto 
the BH occurs relatively frequently in nearby galaxies. 
However, the AGN fraction does not reveal how quickly 
the BHs are growing, which requires a measurement of 
the mass accretion rate and Eddington ratio (see Eqn. 1 
& 6 and Fig. 1). Using a broad range of methods to 
estimate BH masses and adopting the [O iii]/l5007 or 
[O iv]/125.9 i-im luminosity as a proxy for the mass- 
accretion rate, the majority of AGNs in the local Uni- 
verse have Eddington ratios of /I » lO^^'-lO"^, with a 
small tail towards higher Eddington ratios (see Fig. 9 of 
.Ho.2008t [Goulding et alTpOlOj see also Fig. 1). Since 
the growth time (i.e., the time to double in mass) of a BH 
with an Eddington ratio of A ^ 10"^ is comparable to the 
age of the Universe (see Eqn. 10 in §4.3), these range of 
Eddington ratios are modest and suggest that the major- 
ity of the massive BHs in the local Universe must have 
grown more rapidly in the past. However, these results 
can also be a bit misleading since there is a small but 
very significant tail of AGNs with much higher Edding- 
ton ratios. Indeed, x 50% of the mass accretion onto 
BHs in the local Universe is found to occurred at high 
Eddington ratios of A > 0.1, despite these systems only 
comprising ^ 0.2% of the optically identified AGN pop- 
ulation ([Heckman^etaE] |20^| ! These results therefore 



^Radio-loud AGNs are typically identified on the basis of lumi- 
nous radio emission (e.g., Li.4ghz i lO""* W Hz"' ) which is not due to 
star-formation activity. By comparison, the term "radio-quiet AGNs" 
refers to the majority of the X-ray, optical, and IR selected AGN pop- 
ulations, which are not luminous at radio wavelengths. 
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suggest that, while the vast majority of the massive BHs 
in the local Universe are growing slowly (potentially in 
the optically thin accretion mode), at least half of the 
mass accretion in the local Universe occurs is the opti- 
cally thick accretion-disc mode (see §2.3 and Fig. 3). 



AGN activity is transient and therefore the measured 
Eddington ratio of an individual object only provides a 
snapshot of its growth path. We clearly cannot watch 
an individual BH grow over Myrs of time to deter- 
mine its time-average growth but we can measure the 
volume-average Eddington ratio of a group of systems, 
which will provide an average (and therefore typical) 
BH growth rate. On the basis of optical spectroscopy 
from the SDSS, |Heckman et al.1 ( |2004 ) showed that the 
volume-average Eddington ratio for optically identified 
AGNs in the local Universe is inversely proportional 
to the BH mass; see Fig. 1 for a graphical representa- 
tion of this result. The volume-average BH growth time 
for the lowest-mass BHs (Mbh ~ 10^ Mq) is compara- 
ble to the age of the Universe, implying that their BHs 
have been growing steadily over cosmic time; |Gould-| 
ing et al. ( 2010| l extended this work to lower BH masses 
and showed that the most rapidly growing BHs in the lo- 
cal Universe reside in late-type galaxies with low-mass 
BHs (Mbh ~ 10''-10^ Mq), many of which lack the 
clear optical signatures of AGN activity. By compari- 
son, more massive BHs with Mbh i 10^ M© are typi- 
cally growing x 1-3 orders of magnitude more slowly 
than the lower-mass BHs and must have undergone their 



dominant growth phases at higher redshifts. Kauffmann 



& Heckman ( 2009 1 extended these results by providing 
evidence for distinct regimes of BH fuelling in the local 
Universe. They showed that optically identified AGNs 
with significant star formation in the central kpc have 
a lognormal distribution of Eddington ratios (with the 
peak at an Eddington ratio of /I ~ 10"^) while more 
quiescent AGNs were found to have a power-law distri- 
bution of Eddington ratios, where the BH accretion rate 
depends on the stellar spheroid mass and the ages of 
the stars in the bulge. The active and quiescent systems 
are broadly split on the basis of BH mass: the active 
systems have Mbh ~ 10*^ Mq while the quiescent sys- 
tems have Mbh ~ 10^ Mq. On the basis of these results. 



Kauffinann & Heckman ( 2009 1 argued that the BHs in 



the active systems are driven by a plentiful supply of 
cold gas while the BHs in the more massive quiescent 
systems are driven by stellar-mass loss in evolved stars, 
having presumably consumed their cold-gas supplies at 
higher redshift. 



3.2.2. Is black hole and galaxy growth connected? 

The BH-spheroid mass relationship in the local Uni- 
verse suggests that AGN activity should be found in 
spheroid-hosting galaxies. Indeed, the majority of opti- 
cally identified AGNs reside in comparatively early type 



galaxies (E-Sb galaxy morphologies; e.g.. Ho 2008 
[Gadotti & Kauffmann||2009] ). However, a substantial 
fraction of the AGN population is also hosted in late- 
type galaxies (Sc-Sm). This was a somewhat unex- 
pected result since late-type galaxies host pseudo bulges 



rather than classical spheroids (e.g., Kormendy & Ken- 



nicutt 



2004| l, and are therefore not expected to typically 



host massive BHs (since the BH-spheroid mass rela- 
tionship is defined for galaxies with classical spheroids). 
Focused investigations have indeed shown that there 
are no significant correlations between BH mass and 



pseudo bulge luminosity (e.g., [Greene et al. 2008 Jiang 



et al. 201 1 ; Kormendy et al. |201 1 1. The identification 



of growing BHs in late-type galaxies therefore suggests 
a different growth path to systems that host classical 
spheroids; for example, it is possible that these late-type 
systems have never undergone a galaxy major merger, 
and a non-negligible fraction of the BH mass may be 
from the massive BH seed formation (e.g., Kor mendy &] 
Kenmcutt|[2004t [Volonteri & Natarajan„2009, see §4.3 



for the discussion of the formation of massive BHs). 
AGNs hosted in late-type galaxies are often not revealed 
at optical wavelengths and mid-lR spectroscopy or hard 
X-ray observations are required to unambiguously iden- 
tify them (e.g., "Ho lOOS"; Satyapal et al.' "^OOSl [Gould- 



ing & Alexander, ,2009; Koss et al., 2011) : the presence 



of dust and star formation in the host galaxy often dilute 
or extinguish the optical AGN signatures. 

The host-galaxy properties of radio-loud AGNs are 
connected to their optical spectral properties. Radio- 
loud AGNs with low-excitation optical spectra (i.e., 
LINER classifications) are mostly found in massive red 
early type galaxies, while radio-loud AGNs with high- 
excitation optical spectra (i.e., Seyfert and quasar clas- 
sifications) are typically hosted in less massive gas-rich 



galaxies (e.g., Willott et al. 2001 Chiaberge et al. 



|2005||Evansetal.| 20061 IHardcastleetal. 2006 2007 


Baldi & Capetti"2008'i'Kauff^mann et al.' 


2008';'Smolcic 


.2009; Smolcic et al.,,2009, Herbert et al. 


2010,, Smol5i6l 



& Riechers 201 li. To first order the high-excitation 
radio-loud AGNs are high-accretion rate systems and 
are believed to evole strongly with redshift, a behaviour 
that would be similar to radio-quiet Seyfert galaxies 
and quasars; see §3.3.1. By constrast, low-excitation 
radio-loud AGNs are low-accretion rate systems, evolve 
slowly with redshift, and are the dominant radio-galaxy 
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population in the local Universe (e.g., [Sadler et al. 



2007 Smolcic et al. 2009 1. The low-excitation sys 



tems are the dominant radio-loud AGN population at 
modest radio luminosities (Li.4ghz ~ 3 x 10^^ W Hz ') 
and also dominate the global energetic output from ra- 
dio galaxies (e.g., Hardcastle et al. 2007| l; these systems 
are typically classified as FR I radio galaxies ( Fanaroffj 
|& Rileyl |1974[ ) and are important for the "radio-mode" 
AGN feedback (e.g., [Merloni & Heinz|2008; see §5.1). 
The high-excitation systems are typically classified as 
FR II radio galaxies and typically have the highest radio 
luminosities but are comparatively rare at all redshifts. 

The tightness of the BH-spheroid mass relationship 
suggests that galaxies hosting classical spheroids grew 
in concert with their BHs (i.e., for every ^ 1000 units of 
star formation there is ^ 1-2 units of BH accretion). As 
shown in §2.2, a direct link between star formation and 
AGN activity on a; 10 pc scales may be expected given 
that both processes are driven by the availability of cold 
gas on nuclear scales. However, a connection between 
AGN activity and star formation on the kpc scales of 
the galaxy spheroid is not expected a priori since the 
vast difference in size scale would preclude a direct 
causal link. Indeed, although the host galaxies of many 
AGNs are undergoing star-formation activity, there is 
large scatter in the observed AGN: star formation ratio 



(e.g., [Netzer etaLl [2007) [Wild et al.[[2007j Baum et all 



[2010(|Diamond-Stanic & Rieke|[20ri| i7Tlie scatter de- 
creases when the star-formation rate is measured over 
< 1 kpc scales (e.g., Wild et al. 2007 [ Diamond-Stanic 



& Rieke 2011 1, and is predicted to decrease yet fur- 
ther at smaller spatial scales, where a causal connection 
is expected between the gas that drives star formation 



and the gas that fuels the BH (e.g., Hopkins & Quataert 
[2010| l. Despite this, using optical spectroscopy from the 
SDSS, [Heckman et al. (2004) showed that the volume- 
average galaxy-BH growth rate is broadly consistent 
with that expected from the BH-spheroid mass relations 
(a factor of ~ 1000) across kpc scales (i.e., the typical 
size scales corresponding to the SDSS fibres), for a wide 



range of BH masses (see Fig. 5 of Heckman et al. 2004 1. 

The fraction of galaxies hosting AGN activity is 
correlated to the IR luminosity (8-1000 yum) or star- 
formation rate, indicating a connection between the duty 
cycle of BH growth and the growth of the galaxy. On 
the basis of optical spectroscopy alone, the AGN frac- 
tion rises from a; 5% for galaxies with Ljr < 10' ' Lq 
to ~ 15% for Luminous IR Galaxies (LIRGs; Lir a 
10"-10'2 Lo) and a 25% for Ultralu minous IR Galax- 
ies (ULIRGs; Lir > lO'^ Lq; e.g.. 



Kim et al. 



1998 



[Veilleux et ai:][T999l [Tran et al.[[2001| l. However, op 
tical spectroscopy does not identify all of the AGNs and 



these AGN fractions increase to a 50-80% in the LIRG 
and ULIRG populations when including mid-IR spec- 
troscopy and X-ray data for the identification of AGNs 



(e.g.. 


Alexander et al. 2008a| Lehmer et al. 


2010 Nar- 


dini et al.| 20101 Alonso-Herrero et al. 20 


1 1 1 INardini 


Risaliti 201 1 1. These results show that the BH grows 



almost continously during periods of intense star forma- 
tion. 

3.2.3. How is AGN activity triggered? 

Central to our understanding of the growth of BHs 
is the determination of the large-scale physical mecha- 
nisms that trigger the gas inflow towards the BH - either 
external (galaxy mergers or interactions) or internal (gas 
instabilties, galaxy bars, etc; see §2.1 and Fig. 2) pro- 
cesses. Internal processes and galaxy interactions are 
often referred to as "secular evolution" and can be fur- 
ther divided into internal secular evolution and external 
secular evolution, respectively (see Fig. 1 of [Kormendy[ 
[& Kennicutt[[2004] l. An accurate determination of the 
fraction of AGN activity driven by these different pro- 
cesses is non trivial since it depends on the depth of the 
data, the spatial resolution of the data, the assumptions 
used on how the observed data relates to the trigger- 
ing mechanism, and how the AGN activity is identified. 
Furthermore, given the gas-inflow times, there can be a 
significant delay between the initial triggering event and 
the onset of AGN activity (=s 50-500 Myrs; e.g.,|Davies 



[eraLl[2007i fSchawinski et"ai[2009i [WdeFal 



20101; 



see §2.1-2.2. These compUcations mean that it can be 
difficult reconciling the results from one study with the 
results form another study, let alone accurately measur- 
ing the true fraction of AGN activity that is triggered by 
external or internal processes. 

For example, even in the local Universe there are 
large differences between published studies. Using mor- 
phological classifications of SDSS galaxies from the 
Galaxy Zoo project (Lintott et al. 201 1\ , Darg et al. 
( 2010| l estimated that only ^ 1% of AGNs reside in sys- 



tems clearly undergoing a major merger. By compari- 
son, on the basis of a hard X-ray selected AGN sample 
(from the Swift-BAT survey; iTueUer et al.|2008l[20roi l. 



Koss et al. (2010 1 estimated the rate of major mergers 
to be 20%; the fraction rises to a 50% if galaxy in- 
teractions are also included. Part of the discrepancy in 
the estimated rate of major mergers between these stud- 
ies may be due to the AGN identification method, since 
the Swift-BAI AGN survey finds AGNs that lack un- 



ambiguous optical AGN signatures (e.g.. Winter et al 



[2009 ) and may probe higher intrinsic luminosities than 
the optical samples ( [Koss et "aL] [201 1) . Still, a fac- 
tor a 20 difference seems difficult to reconcile from 
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the AGN selection method alone. In another example, et al. 201 l| l, and are therefore unlikely to have ever 



Liu et al. (201 1 1 looked at the frequency of AGN pairs 
within ^ 5-100 kpc in the SDSS and found that ^ 4% 
of optically identified AGNs are found in pairs. How- 
ever, only X 30% of these AGN pairs clearly showed the 
morphological features of a major merger, even though 
the host galaxies appear to be in the process of a galaxy 
merger or interaction. The depth of the imaging data 
used to identify the morphological signatures of major 
mergers can clearly have a large effect on the results. 



Using deep imaging data, Schawinski et al. (2010 1 iden- 
tified the faint morphological signatures of galaxy ma- 
jor mergers in > 50% of blue early type galaxies (these 
weak features often referred to as "shells"; |Malin &] 



Carter] 1983[|Hernquist & Spergel|1992}|Turnbull etH 
1999[ l, which were not apparent in shallower data. 

Spectroscopy can provide an alternative route to esti- 
mating the fraction of systems in major mergers by iden- 
tifying systems with evidence for close-pair AGN activ- 
ity. The identification of double peaked emission lines 
(broad or narrow) with significant velocity offset can in- 
dicate a binary system or a BH gravitational recoil (e.g., 
Volonteri '200T; Komossa et al. '2008'; Colpi & Dotti' 
2009, .Dotti et al., .2009, .Civano et al., .2010, .Robinsoni 



et al.||2010i[Smith et al.||2010l|Rosario et aLl[20ni[Barn 



rows et air '201 1'); however, see [Dotti & Ruszkowski| 
( 2010[ l and ^Shen et al.,(2011 ) for alternative interpreta- 
tions of the presence of double-peaked emission lines. 
On the basis of the incidence of double-peaked emis- 
sion lines, Comerford et al.| ( |2009| l estimated a ^ 30% 
merger fraction, similar to the upper end of merger frac- 
tions estimated from imaging data. 

These suite of studies highlight the difficulties in ac- 
curately quantifying the triggering mechanism of AGN 
activity. Major mergers clearly occur and trigger AGN 
activity but there are significant uncertainties in the in- 
terpretation of the data; indeed, systems with classi- 
cal spheroids may have undergone at least one major- 
merger event at some time during the past (e.g., Bour- 



|naud et al.| [2005] [Hopkins et al.| |20T0l l. Conversely, 
AGN activity is also clearly often triggered by galaxy 



interactions and internal processes (e.g., [Malkan et al. 



[T9981 [Kuo et"aLl [2008) plison et al.[ [201 l| l. An exam- 
ple of an AGN population where the majority of the BH 
growth may be entirely driven by secular processes are 
narrow-line Seyfert 1 galaxies (NLSls), a subset of the 
Seyfert 1 population with comparatively narrow broad 
emission lines (p ermitted broad-line widths o f ~ 500- 
3000kms-i- 



e.g. 



Osterbrock & Pogge 



1985 I. The ma- 



jority of NLSls reside in galaxies with pseudo bulges, 
in contrast to the classical spheroids hosted by typical 



undergone a major merger (e.g., Kormendy & Kenni 
[cutt[[2004) . Pseudo bulges have more angular momen 
tum than classical spheroids, which inhibits the gas ac- 
cretion on small scales. It is therefore significant that 
NLSls are more likely to reside in barred spiral galax- 
ies than typical Seyfert galaxies (Crenshaw et al. 2003| l 
- without the galaxy bars, the host-galaxy gas is un- 
likely to be efficiently driven into the central kpc region 
of NLSls (contrast with the lack of a strong correlation 
in the overall population; see §2.1). 

3.3. The distant growth of massive black holes 

3.3.1. Evolution of the AGN population 

Extensive high-quality observations have provided 
valuable insight into the processes of BH growth in 
the local Universe. However, the vast majority of the 
mass accretion onto BHs occurred at higher redshift 
(~ 95-99% of the integrated BH growth has occurred at 
z > .1; e.g.,|Marconi et aL|2004, Shankar et al..2004l 



2009} [Hopkins et al.|2007[|]\4erloni & Heinz|2008[[A]rd 



et al.[|20l()| l, and we must therefore look to the distant 
Universe to understand when, where, and how today's 
massive BHs grew. Optical quasar surveys established 
more than four decades ago that luminous AGNs were 
orders of magnitude more common at z Z 1-2 than 
z ^ (e.g., [Schmidt[ [19681 [Schmidt & Gre"enl [T983 



'Hartwick & Schade[[1990[ [Boyle et al.[ [2000| [Rich^ds| 
,et al. ,2006) . However, the optical photometric selection 
method used to identify quasars is only effective in se- 
lecting unobscured AGNs that are bright enough to out- 
shine the host galaxy. Optical quasar surveys therefore 
provide limited constraints on the properties and evolu- 
tion of the majority of the AGN population; i.e., either 
obscured AGNs or lower-luminosity AGNs. 

Currently the most efficient and near-complete selec- 
tion for distant AGNs is made using X-ray observations. 
X-ray observations provide an almost obscuration- 
independent selection of AGNs and, since star forma- 
tion activity is comparatively weak at X-ray energies, 
can identify even low-luminosity systems. For example, 
the deepest X-ray surveys (the Chandra Deep Fields; 
CDFs; [Brandt et al.'200r; Giacconi et al.'2002 ; Alexan-[ 
[der et al.t2003b ; Luo et al. 2008 , Xue et al. 2011) can 
detect low-luminosity AGNs similar to NGC 405 1 (see 
Fig. 1 and also Fig. 6 of Brandt & Hasinger"2005) out 
to z ^ 1-2, even in the presence of large amounts of ab- 



sorption (up-to A^h ~ 10^^-10^"* cm ^; e.g.. 


Tozzi et al. 


'2006'; 'Raimundo et al.'2010 ; Alexander et al.|2011l Co- 


mastri et al. 2011 


Feruglio et al. 20lTll; AGNs ^ 1-2 



Seyfert galaxies (Mathur et al. 201 1[ Orban de Xivry orders of magnitude brighter can be identified to z > 6, 
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provided a sufficient number of objects exist in the com- 
paratively small survey volumes. Often the significant 
challenge in the identification of distant X-ray selected 
AGNs is an accurate measurement of source redshifts 
(with spectroscopic or photometric data; e.g., Barger" 



et al.||2003b[ [Szokoly et al.|[2004t [Zheng et al.,|2004t 
Cardamone et al.||2010[ |Luo et aLpOlOtlSalvatoeraL] 



201 l| l, since the optical/near-IR counterparts for many 



of the AGNs are very f aint (e.g., [Alexander et al. 2001 
[Mainieri et aL] |2005b| [Rovilos et al.| |2010| l. Further- 
more, even the deepest X-ray surveys miss the most 
heavily obscured luminous AGNs where even hard X- 
ray photons are absorbed; selection techniques using 
optical spectroscopy, IR, and radio data are starting to 
identify large numbers of these systems (e.g., Donley 



let al.||2005t|baddi et al.i[2007al [Alexander et al.||2008b] 
[Fiore et al.| [20081 [Hickox e taT '20091 [Yan et al.| [2011] 
[Juneau et al.| [201 1| [Luo et al... 20lT). However, de- 
spite these challenges, from a combination of X-ray ob- 
servations across a broad range of the flux-solid angle 
plane (e.g., see Fig. 1 of Brandt & Hasinger|2005 1 



are starting to piece together a more complete picture of 
the evolution of AGN activity across cosmic time. 

The evolution in the space density of high-luminosity 
AGNs broadly tracks that found from optical quasar 
surveys: the space density of luminous AGNs with 
i2-i0kev > 10'^'*-10'*^ erg s"' peaks at Zpeak = 1-2 
and drops by ^ 2-3 orders of magnitude to z =s and 
=s 0-1 order of magnitude to z * 5 (e.g., Fiore et al.' 



20031 [Ueda et al. 


[20031 [Shankar et al.| ,2004, .2009, 


Barger et al. 


2005 


Hasinger et al. 2005 [ Hopkins etal. 


20071 Silverman et al. 2008a| Brusa et al.| 2009a| Aird 


etal. 2010) 

1 


The constraints are currently very uncer- 



1 1 1 1 

tain for z > 5 AGNs (e.g.. 


Fan et al. 


2001 


Barger et al. 


2003a Willott et al. 2010)1; see §4.3. By contrast, lower 



luminosity AGNs evolve more slowly and their space 
density peaks at lower redshifts than high-luminosity 
AGNs: for example, the space density of moderate- 
luminosity AGNs with L2-iokeV ~ lO'^-'-lO'*'* erg s"^ 
peaks a.tz~ 0.5-1.0 and drops by x 1-2 orders of mag- 
nitude to z ^ and z ~ 5, respectively. This differen- 
tial redshift evolution of the AGN population is com- 
monly modelled as luminosity dependent density evo- 
lution. The same behaviour is also found for optically 



selected and radio selected AGNs (e.g., J^chards et al. 



et al., .2009j n 

Obscured AGNs are found to trace the 



2005[ [B ongiorn o et al.[ [20(j7[ [Hopkins et al.[ [2007 
Rigby et al. i2008| [201 11 [Croom et al.| [20091 [Smolcic 



same evolution as unobscured AGNs of the same lu- 



'Analogous results are also found for the star-forming galaxy 
population, where the space density of high-luminosuty star-forming 



minosity, although there is tentative evidence for an in- 
crease in the obscured AGN fraction with redshift (e.g.. 



La Franca et al.[|2005| 


Ballantyne et al.[[2006t [Treister 


& Urry[[2006| Hasinger 


|2008 


however, see also Aky- 


las et al.|2006 and Dwe 


ly & 


'age|2006 1, which would 



be expected if the nuclear regions of distant AGNs are 
more gas rich than lower-redshift AGNs. 

The integrated growth of BHs is dominated by sys- 
tems around the knee of the AGN luminosity func- 



tion (e.g., Hopkins et al. 2007 1, which for X-ray de- 
tected AGNs is L2-iokeV * 10^ erg s^': ^ 75% of the 
growth of BHs has occurred in luminous systems with 
La-\QkeV ~ lO'^-'-lO'*^ erg s"' (i.e., straddling the tradi- 
tional Seyfert galaxy/quasar threshold; see footnote 1 of 



Alexander et al. .2008b ). Approximately 30-50% of the 



integrated growth of BHs has occurred at comparatively 
low redshifts of z < 1, ~ 35-45% has occurred at z ~ 1- 



2, and ^ 15-25% at z > 2 (e.g., [Marconi et al. 


2004 


Shankar et al. 


2004 2009 Hopkins etal. 2006a| 


2007 


Silverman et a 


I. 2008a Airdetal. 2010 1. Optically se- 



lected quasars account for ^ 35-50% of the integrated 
BH growth and radio-loud AGNs account for t 10% 
of the integrated BH growth. However, while compris- 
ing the minority of the overall AGN populations, both 
of these sub populations are important for the growth 
of BHs and galaxies: optically selected quasars appear 
to represent a rapid growth phase of massive BHs (see 
§4.2) and radio-selected AGNs appear to have played a 
significant role in the formation and evolution of galax- 
ies (see §5.1). 

The origin of the differential redshift evolution of the 
AGN population is presumably due to a decrease in the 
overall cold-gas fuel supply in the nuclear regions of 
galaxies, at least for the lower-redshift drop off in space 
density. The global drop off in space density at z > 2 
may be limited by the maximum BH accretion rate for 
comparatively low-mass BHs (i.e., the AGN luminos- 
ity cannot be higher than the BH Eddington limit; e.g., 
Merloni & Heinz 2008^ see Eqn. 6). However, it isn't 



immediately clear whether the decrease in the fuel sup- 
ply is across all systems or only for systems in specific 
environments or BH mass ranges. Perhaps the strongest 
discriminator between these different scenarios comes 
from studies of nearby AGNs, where the long growth 
times of the most massive BHs at z ~ 0.1 (Mbh > 
10^ Mq) implies that they must have grown more rapidly 



galaxies peaks at higher redshifts than lower-luminosity systems (e.g., 
|Le Floc'h et al.| |20()5 Perez-Gonzalez e t al.|[2005) . This behaviour 
is typically referred to as galaxy "downsizing" since there is clear ev- 
idence for a decrease in the mass of luminous star-forming galaxies 
with decreasing redshi ft (e.g.,|Cowie et al.||1996[ [Ballantyne et al.| 
[2006|[Bell et"ar[[2005[|Juneau et al.[[2005[[Bundy et al.[[2006| . 
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at higher redshifts (e.g., Heckman et al. 2004] l; see 
§3.2.1. Direct measurements of the Eddington ratios of 
the BHs in z < 2 quasars (BH masses estimated using 
the virial estimator; see §4. 1 ) provide further support for 
this view by showing a decrease in both the character- 
istic active BH mass and Eddington ratio with decreas- 



ing redshift (e.g., McLure & Dunlop J2004 Greene & 
pol[2007l|Netzer & Trakhtenbrot||2007| l. Lastly, on the 
basis of indirect BH mass estimates (based on the host- 
galaxy luminosities or velocity dispersions and the local 



BH-spheroid mass relationships; e.g., Magorrian et al. 
[T9981 [Gebiiardt et al.|2000, .Ferrarese & Merritt..2000, 
ITremaine et"aL][2002l [Marconi et aTHIOO?) !, the BHs 
of z ^ 0.3-1.5 X-ray AGNs appear to be massive and 
are growing more rapidly than similarly massive BH 
in the local Universe (Mbh ~ 10^-10^ M q and typical 
Eddington ratios of A ^ 10^^-10 '; e.g., Babic et al. 



2007 Ballo et al. 2007 Rovilos & Georgantopoulos, 



20071 lAlonso-Herrero et al.||20081 [Hickox et al.||2009r 
Raimundo et al.|2 010; Simmo ns et al.|2011| l. 

The current observational data therefore suggest 
"downsizing" in the active BH mass with decreasing 
redshift, in general agreement with that found for the 
galaxy population (see footnote 7). However, the cur- 
rent observations provide comparatively limited con- 
straints in isolation since accurate determinations of the 
BH mass and Eddington ratio are only available for a 
small fraction of systems (i.e., AGNs with broad emis- 
sion lines; see §4.1), and the constraints are very lim- 
ited for z > 1.5 AGNs, where the most massive BHs 
are predicted to have been most active. To gain more 
detailed insight we can appeal to models and simu- 
lations, which provide solutions to the growth of the 
BH population from a broad suite of observational and 
physical constraints (e.g., BH mass density and dis- 
tribution, evolving AGN luminosity functions. X-ray 
background spectrum; e.g., Marconi et al.||2004[ Mer- 
loni' 2004; Shankar et al. 2004 2009' 'Hopkins et al 
2006a, 2008 , Malbon et al. 2007; Di Matteo et al. 2008 



Somerville et al.|2008||Fanidakis et al.|201 l|pOT2] i. De- 
spite a broad range of different approaches, from sim- 
ple analytical/synthetic models to more detailed semi- 
analytical models and N-body/hydrodynamical simula- 
tions some clear trends have emerged from these stud- 
ies; (1) the average Eddington ratio for BHs decreases 
with redshift, from rapidly growing systems atz~ 2-3 
to comparatively slow growing systems at z < 1, and 
(2) the characteristic "active" BH mass decreases with 
redshift. However, there are significant differences be- 
tween the studies in terms of the build-up of the BH 
mass function with redshift. For example, in the semi- 



mass end of the BH mass function (Mq i 10^ Mq) is 
essentially in place by z ~ 2 while the high-mass end is 
predominantly built by BH mergers of lower-mass BHs 
and modest mass accretion rates through to the present 
day. By comparison, in the AGN synthesis models of 
[Marconi et aT] ( |2004l l and [Merloni & Heinz[ ( |2008l l the 
lower-mass end of the BH mass function is only in place 
by z ~ 0.1-0.6. Of course there are large differences in 
the complexities, scope, and assumptions made between 
these different models but they serve to illustrate some 
of the uncertainties that remain in modelling the evolu- 
tion of the BH mass function. 

Radio-loud AGNs are the minority AGN population 
at all redshifts and the majority of the growth of BHs has 



occurred in radio-quiet AGNs (e.g., Merloni & Heinz 



2008 Cattaneo et al. 2009 1. However, a large frac- 



tion of the radio emission associated with AGN activ- 
ity is produced by jets and lobes, which are powered 
by synchroton emission and shocks and therefore large 
amounts of kinetic/mechanical energy. If this kinetic 
energy is able to couple effectively to the gas in the host 
galaxy or larger-scale environment then it can prevent 
the gas from cooling and forming stars; see §5.1. The 
evolution of the radio luminosity density with redshift 
therefore provides direct constraints on the volume- 
average heating rate from AGN activity (e.g., [Croton 



leTaTI [20061 [Heinz et al.| |2007i [Lehmer et al. 



2007^ 



Merloni & Hemzj [20081 ICattaneo et aLj [2009j [SmoISigl 
let al.| |2009| [La Franca eTall |20T0jl. The conversion 
from radio luminosity to mechanical power is uncertain 
and relies on converting the synchrotron jet luminosity 



into a kinetic energy (e.g., Dunn & Fabian 2004 ; ^est| 



et al. 2006 [Rafferty et al ,2006, .Heinz et al., ,2007 



Birzan et al.] [2008| [Cavagnolo et alTl [2010| l. However, 



test] 
07) 



on the basis of the current conversion factors, the over- 
all heating rate from AGN activity is predicted to be 
broadly flat over z ~ 0^, with low excitation AGNs 
dominating the overall heating, particularly at z < 1 
when radio-loud high-excitation AGNs are rare (e.g., 
[Merloni & Hemz] [20081 [Cattaneo et al.[|2009l l. There is 
evidence for a sharp drop in the heating rate at z < 0.5, 
which would indicate a weakening role of AGN activity 



towards the present day (e.g., Kording et al. 2008 La 



Franca et al. 2010 1. However, there are also significant 



analytical model of Fanidakis et al. (2012l, the lower- 



differences between studies in the estimated evolution 
of the AGN heating rate, which are due in part to how 
the current radio luminosity functions are extrapolated 
down to low luminosities. 

3.3.2. The host galaxies of distant AGNs 

The strong evolution in the AGN population with red- 
shift illustrates that conditions in the distant Universe 
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were different to that seen locally. How do the host 
galaxies of distant AGNs compare to the host galaxies 
of AGNs in the local Universe? Accurate measurements 
of the host-galaxy properties of distant AGNs can be de- 
termined from high-spatial resolution imaging, where 
the AGN component can be modelled as a point source, 
or from photometry for obscured and low-luminosity 
systems when the host-galaxy emission dominates at 
rest-frame optical-near-lR wavelengths. 

The concensus view is that the majority of distant 
AGNs are hosted in massive galaxies. The host galax- 

1 0"^^ erg s~ ^ at z < 3 are 



et al. 



(0.3- 3) xlO" Me 
2007; Ballo et al. 



(e.g.. 


Akiyama 2005 


Babic 


2007; Alonso-Herrero et al.| 



2008 , Bundy et al., ,2008, ,Brusa et al., ,2009b, ^icko5c| 



et 



aL^ ' 12009 



Xue et al. 



2010[ l. Lower-mass systems 



(M, * (0.3-2) xlO"' Mo) are also identified but appear 
to comprise the minority population of X-ray AGNs 
(e.g., |Shi et al.| [20081 |Xue et al.| |20T0}; however, se- 
lection and sensitivity effects means that it is often chal- 
lenging to identify distant AGNs in low-mass galaxies 
(i.e., for a given luminosity threshold, a low-mass BH 
needs to be growing at a higher Eddington ratio to be 
identified than a high-mass BH; see Fig. 1). Radio- 
loud AGNs are typically hosted in the most massive 
elliptical galaxies with M» ^ 10"-10'^ Mq: systems 
with low-excitation optical spectra appear to be hosted 
in higher-mass galaxies than high-excitation radio-loud 
AGNs, as expected if the low-excitation systems reside 
in the most massive dark-matter halos and underwent 



their major growth phases at higher redshifts (e.g., Tasse 



etai:i|2008|perbert et al.[[20T0| [Floyd etalj [20101 ); see 
§3.4. By comparison, optically selected quasars are pre- 
dominantly hosted in less-massive elliptical galaxies, 
with a non-negligible fraction found in galaxies with a 



significant disc component (e.g., Dunlop et al. 2003 



Royd et al. 2004 Kotilainen et al. 2007 Bennertet al 
2008', 'Schramm et al. 2008; Tasse et al.[ |2008t |Kod^ 



lainen et al., 2009, .Veilleux et al., 2009| [Floyd et al. 

2oTo| 



Various analyses show that the host-galaxy proper- 
ties of distant AGNs are also similar to those of distant 
massive galaxies. In terms of the Colour-Magnitude 



Diagram (CMB; e.g., |Strateva et aL[[2001[ ), distant X- 
ray AGNs predominantly lie in the "green valley" be- 



tween the "red sequence" and "blue cloud" (e.g.. Nan 



dra et al. 2007; Rovilos & Georgantopoulos , '2007! Sil 



verman et al., 2008b, Hickox et al., 2009, Xue et al. 



2010| [Cardamone et al.| [2010| l, consistent with that of 
inactive galaxies over a similar mass range; however, we 
note that some AGN host galaxies will be red due to dust 



201 1 1. The host galaxy morphologies of X-ray detected 
AGNs are also comparable to those of similarly mas- 



sive inactive galaxies (e.g., Sanchez et al. 


2004 Grogin 


[et al.||2005||Pierce et al.||2007l Gabor et a 


1. 120091 Geor- 


'gakakis et al. 2009 Kocevski et al. 2011 Schawinski 


iet al. 


,201 1|, with a mix of bulge-dominated and disc- 



dominated systems. By comparison, the host galaxies 
of distant radio-loud AGNs typically lie in the "red se- 



quence" of the CMD (e.g., |Tasse et al. 2008} Hickox 



et al. 2009 1, as expected for massive elliptical galax- 
ies. The morphologies and structural properties of dis- 
tant radio-loud AGN are also typically consistent with 
those of distant elhptical galaxies (e.g., Dunlop et al. 
"^OOT; M cLure et aL] [20041 [Floyd et al.| |2010| [Herbert 
|etal.,,201lT [ 

The similarity in the host-galaxy properties of dis- 
tant AGNs and coeval inactive galaxies suggests that the 
AGN fraction provides an estimate of the BH-growth 
duty cycle. The fraction of massive galaxies hosting X- 
ray AGN activity out to z 2-3 i s ^ 10-20% for sys- 
tems with Lx i 10"*^ erg s"^ (e.g., [Bundy et al.[ [2008 



Xue et al.[ [2()T(j| ). The AGN fraction drops to 
10« 



1-5% 

with some evidence 



for AGNs with Lx ~ iU"" erg s 
for an increase with redshift (e.g., [Xue et aLl [2010[ 
Georgakakis et al. 201 l| l, comparable to that found in 



the local Universe but at higher BH and stellar masses 
(see Fig. 4). These results indicate that BH growth is re- 
current and has a long duty cycle. The fraction of galax- 
ies hosting radio-loud AGN activity out to z « 1.3 is 
also comparable to or slightly higher than that found for 
nearby galaxies, indicating that the duty cycle of radio- 
loud AGN activity has been relatively constant over at 
least the past ~ 9 Gyrs (e.g., Tasse et al. 2008,,Smolcic[ 



[etaLl[2009l l. 

How is the distant AGN activity triggered? Contrary 
to initial expectations, only a comparatively small frac- 
tion of the X-ray AGNs out to z ^ 3 clearly reside in 



galaxy major mergers (s; 15-20%; e.g.. 


Sanchez et al. 


2004, Pierce et al. 2007| Georgakakis et a 


. 2009. Cister- 


nas et al.|201 IJ Kocevski et al.|201 1 comparable to that 



found for X-ray selected AGNs in the local Universe 
(e.g., Koss et aL][2010 see §3.2.3). However, the host 



galaxies of up-to ^ 45% of the X-ray AGNs show some 
evidence for disturbed morphologies, suggesting that a 
considerable fraction of the distant AGN activity could 
be driven by external processes (i.e., either galaxy ma- 
jor mergers or interactions; e.g., Kocevski et al.[[20TT 



obscuration (e.g., Cardamone et al. 2010 Rovilos et al. 



[Silverman et al.[[2011} . As noted in §3.2.3, there are 
potential complications in the interpretation of the host- 
galaxy morphologies, making it difficult to determine 
the true fraction of distant AGN activity driven by ex- 
ternal and internal processes. It is therefore important 
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to note that, irregardless of these compHcations, there 
are no clear differences in the morphological proper- 
ties of AGNs and comparably massive inactive galaxies, 
suggesting that the mechanisms that trigger distant BH 



growth also trigger distant galaxy growth (e.g., Cister- 



|nas et al.|[20ni|Kocevski et al.|[20n]|Schawinski etal 
|20TT] i. 

The fraction of optically selected quasars with evi- 
dence for morphological distortions varies from study 
to study but is always large 30-100%) and, im- 
portantly, is found to be higher than comparably mas- 



sive distant galaxies (e.g., BahcaU et al. 



1997 Canahzo 



& Stockton! [200T] [Dunlop et al.| |2003| |Guyon et al. 



2006; Benn ert et al.| |2008[ |Urrutia et aL[ |2008| l. The 
overall fraction of radio-loud AGNs with evidence for 
disturbed morphologies is also higher than that found 



for distant massive galaxies (e.g., de Koff et al. 1996 



Dunlop et aL| [2003| [Rinos Almeida et al.| |2011b.a) . 
However, [Ramos Almeida et al. ( 201 lb|a| i "showed that 



the disturbed morphology fraction of radio-loud AGNs 
depends on the optical spectral properties: > 90% 
of the high-excitation systems have disturbed host- 
galaxy morphologies (i.e., similar to the optically se- 
lected quasars), as compared to only » 30% of the 
low-excitation systems. Overall these results therefore 
suggest that typical X-ray AGNs and low-excitation 
radio-loud AGNs evolve along with the coeval massive 
galaxy population, while quasars (either radio quiet or 
radio loud) are often triggered by external processes and 
evolve more rapidly (see §4.2 for further discussion of 



quasars), as predicted by some models (e.g., Hopkins 
etaL] [2008) [Hopkins & Hernquist|[2009] l. 



3.3.3. The connection with star formation 

There are a number of pieces of evidence that con- 
nects the global evolution of star formation with that 
seen for AGN activity: (1) the differential redshift evo- 
lution of the AGN population is also found for the 
star-forming galaxy population (see footnote 7), (2) the 
redshift distribution of the most strongly star-forming 
galaxies traces that seen in the optical quasar popula- 



tion (e.g., [Chapman et al.| [20(j5||Wardlow et al.||2011[ ), 
and (3) the overall shapes of the volume-average mass 
accretion and star-formation histories are broadly simi- 



lar (e.g., Boyle & Terlevich 1998 


[Merloni et al. 2004 


[Silverman et al. 2008al Aird et al 


.| 20101. Such a con- 



nection is not unexpected given the tightness of the BH- 
spheroid mass relationship in the local Universe and the 
association between AGN activity and star formation in 
nearby galaxies; see §3.2.2. However, given the differ- 
ent conditions between the local and distant Universe, 
how does the connection between distant AGN activity 



and star formation compare to that found in the local 
Universe? 

The volume -average mass accretion and star- 
formation histories are similar up to z a; 1.0-1.5, when 
the mass accretion history is scaled up by a factor 
^ 5000 (e.g., |Silverman et al.[[2008a|[Aird et al.[[20T0l i. 
However, there are significant diff^erences at z > 2, 
where the slope of the mass accretion history is much 
steeper than that found the star-formation history; for 
example, there is potentially « 1-2 orders of magnitude 
more star formation per unit volume than that found for 
the mass accretion at z a; 6. At least part of the varia- 
tion with redshift could be due to incompleteness in the 
AGN selection (e.g., there may be an increasing fraction 
of X-ray undetected heavily obscured AGNs with red- 
shift; see §3.3.1) but it also possible that the connection 
between AGN activity and star formation was different 
at high redshift. 

There is a general correlation between star-formation 
rate (SFR) and mass accretion rate for individual dis- 



tant AGNs (e.g. , jSchweitzer et alj [2006[ [Lutz et al 
[20081 [ShTetan |2009[ [Silverman et al.| |2()09bt |Trichas| 



eTaT] [2009] [Rafferty et all [lOTT) . However, the cor 



relation is broad, particularly for AGNs detected in X- 
ray surveys where the AGN selection is comparatively 
complete (i.e., up-to five orders of magnitude variation 
for individual X-ray AGNs; e.g., see Fig. 14 of ,Rafferty[ 
et al.|20lT) . Globally the star-formation rates of AGNs 



of all classes are found to increase with redshift, in gen- 
eral agreement with the broad trend found for the over- 



aU galaxy population (e.g., Archibald et al. 


200 1| Ser- 


jeant & Hatzimin; 
2011b, Lutz et al. 


loglou, 2009VMullaney et al. 2010 


2010, Shao et al. 


2010, .Xue et al. 


2010 Mainieri et al. 2011 Seymour et al. 201 If. For 



example, the average star-formation rates of moderate- 
luminosity AGNs (i.e., systems with Seyfert galaxy lu- 
minosities of L2-iokeV > 10"^"' erg s"') are a 40 times 
higher at z a 2-3 than found for similarly luminous 
AGNs at < 0.1 (e.g., [MuUaney et al.[ [SoTTbj ). This 



rate of increase is in good agreement with that found 
for comparably massive galaxies over the same redshift 



range (e.g., Daddi et al. 2007b Pannella et al. 2009 



Rodighiero et al. 2010a Elbaz et al. 201 l| l, suggest 
ing that the presence of an AGN does not have a sig- 
nificant influence on the total amount of star formation 



in the host galaxy (e.g., Shao et aL[ |2010[ [Xue et alT 
2010 Mullaney et al. 2011b| l; however, we note that 



the situtation may be different for rapidly growing BHs 
in starburst galaxies, which may follow a different evo- 
lutionary path (see Fig. 6 in §4.2). 

Taken at face value these results therefore indicate a 
significant increase in the amount of star-formation ac- 
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tivity for a fixed mass-accretion rate with redshift (i.e., 
there was x 40 times more star formation per unit of 
AGN activity at z ~ 2-3 than found at z < 0.1; e.g., 
MuUaney et al. 201 lh\ , which appears inconsistent with 



the concordant BH-galaxy growth ratio implied by the 
BH-spheroid mass relationship. These results can be 
reconciled with the relative growth implied by the BH- 
spheroid mass relationship if, for example, the AGN 
duty cycle increases signficantly with redshift (see be- 
low for current constraints) or if the majority of the star 
formation in the distant Universe occurs in galaxy discs 
rather than galaxy spheroids. 

The AGN fraction is found to be a strong function 
of star-formation rate in distant galaxies. For example, 
the fraction of moderate-luminosity AGNs with Lx > 
10'^^ erg s"' increases from » 3-10% for SFRs x 30- 
200 Ms yr"' (i.e., equivalent to LlRGs) to ^ 10-40% for 
SFRs * 100-500 Mo yr ' (i.e., equivalent to ULlRGs; 



e.g.. 


Alexander et al. 2005 


Laird et al. 2010 


Syme- 


onidis et al."2010'; Xue et a 




2010 Georgantopoulos 


et al.,,2011, ,Rafrerty et al. 201 


1 1. The AGN fraction 



increases further if lower-luminosity AGNs are consid- 
ered (see Fig. 12 oFRafferty et a l.|2011| l. The high AGN 
fraction at the highest SFRs indicates almost continous 
BH growth during vigorous star-formation phases, as 
would be expected during rapid growth phases. How- 
ever, the overall AGN fractions are broadly consistent 
with those found for intense star-forming galaxies in 
the local Universe (e.g., the fraction of nearby ULIRGs 
hosting AGNs with Lx > 10"*^ erg s"' is ^ 40%; Alexan- 
der et al.|2008b i, suggesting that the duty cycle of BH 
growth in star-forming galaxies has remained relatively 
constant over the broad redshift range of z ~ 0-3. This 
result is found despite the large increase in SFR for 
AGNs over z ~ 0-3 and the global increase in the cold- 
gas mass fraction in massive galaxies with redshift (e.g., 
IDaddi et al.[ |20T0l [Ticconi et al.| |20T0l |Geach et al.| 
[20Tl] i. 

3.4. The role of environment 

The BH-growth results discusssed so far are predom- 
inantly for AGNs identified in field environments (i.e., 
the typical regions of the Universe). How does the 
growth of BHs differ as a function of large-scale struc- 
ture and environment? Theories of large-scale struc- 
ture formation predict that galaxy growth is acceler- 
ated in regions of high density (e.g.,'Kauffmann 1996 
De Lucia et al. 2006). The relationship between stel- 



lar age and local environment provides support for this 
hypothesis, showing that the most evolved and massive 
spheroids reside in the highest density regions (clusters) 



2009 1. Large-scale environment is therefore also likely 
to have a significant effect on the triggering of AGN ac- 
tivity. 

In the nearby Universe, radio-quiet and radio-loud 
AGNs preferentially reside in different environments. 
On the basis of the two-point correlation function, 
which measures the clustering strength of selected pop- 
ulations, the dark matter halos of nearby radio-quiet and 
radio-loud AGNs ar e inferred to be ^ 10'^ and 
=a 2 X 10^^ Mo (e.g.,|Mandelbaum et al.||2009||Donoso 



et al. 2010 1. The clustering stength of the optical 
AGNs is also consistent with that found for compara- 
bly massive galaxies. By comparison, radio-loud AGNs 
are preferentially found in galaxy groups and poor-to- 
moderate richness galaxy clusters (Best 2004[ l. For ex- 



ample, measuring environment as the local density of 
galaxies in a radius of 2 Mpc, the fraction of galaxies 
hosting radio-quiet AGN activity is ^ 2 times higher in 
low-density regions than in high-density regions ( |Kauff-| 



mann et al. 2004), while the local galaxy density of 
radio-loud AGNs is 2 higher than that found for radio- 
quiet AGNs ( jKauffmannetal] [2008] ). 



Broadly similar results are found for distant AGNs. 
The clustering of radio-quiet AGNs is weaker than 
radio-loud AGNs out to z ~ 1-2, with implied dark- 
matter halo masses of ^ lO'^-lO'^ Mo and ^ (0.3- 



1.0) xlO"* Mo, respectively (e.g., Li et al. 2006 da 


[Angela et al.||2008 Coil et al.||2009 


GilU et al. 


2009 


IHickox et al. 2009, 2011'; Mandelbaum et al. 


2009 


iKrumpe et al. 2010, Fine et al.2011 


however, also see 
Jence that the clus- 


|Bradshaw et al.|201 1 1, with some evic 



tering strength is dependent on the adopted AGN lumi- 
nosity thresholds (e.g., ,Krumpe et al. 2010) ; see Fig. 5. 
X-ray AGNs out to z ~ \ are found to reside in a broad 
range of environments (from field galaxy environments 
to galaxy groups), consistent with those found for galax- 
ies of the same mass, colour, and star-formation rate 



(e.g., Georgakakis et al. 2008 Silverman et al. 2009a 



|Tasse et al._, 201 1) . The X-ray AGN fraction for galaxies 
found in the field and galaxy groups environments are 



also indistinguishable (e.g., Waskett et al. 2005 Silver 



man et al. 2009a) . By comparison, distant radio-loud 
AGNs are found to typically reside in galaxy groups 
and clusters, where the radio source may be triggered by 



weak accretion of the hot group/cluster gas (e.g., Tasse 



at the present day (e.g., Baldry et al. 2004 Smith et al. 



[etaLl [20081 [Bardelli et al.| |20T0l [Smolcic et al.||20TTF 
see Fig. 2. 

Stark differences in the incidence of luminous AGN 
activity is found between distant galaxy clusters and 
distant field-galaxy regions. The moderate-luminosity 
AGN fraction (L2-iokeV > 10"^"' erg s"') in cluster galax- 
ies rises by an order of magnitude over the redshift range 
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Figure 5 : Illustration of the evolution of dark matter halo mass ver- 
sus redshift for AGN populations. Lines showf the evolution of halo 
mass with redshift for individual halos, based on the median growth 
rates from cosmological simulations jpakhouri et al.|[2010^ . High- 
lighted is the region of maximum "quenching", in which halos tran- 
sition from having large reservoirs of cold gas to being dominated 
by virialized hot atmospheres (e.g., Croton 20091. The gray points 
show halo masses of optically-bright quasars derived from autocor- 
relation me asurements iCroom et al. 2005 upright triangles; Myers 



et al. 


2006 


et al. 


2008 



squa res; |Shen et al.|2007| inverted triangles; |da Angelaj 
stars; |Ross et al.|2009| circles). The colored points show 
the halo masses for radio. X-ray, and infrared-selected AGN at z ~ 0.5 
jHickox et al.||2009) . This figure illustrates that massive rapidly grow- 
ing BHs (optical quasars) are found in the most massive dark matter 
halos that have not yet reached the "exhausted" hot halo regime, while 
rapidly growing lower-mass BHs (as traced by mid-IR samples) are 
found in somewhat lower-mass halos typical of star-forming galaxies. 
All halo masses aie calculated following |Hickox et al.|j20TT} assum- 
ing a cosmology with (£!„,, £1a. fg) = (0.3, 0.7, 0.84). 



0.05 < z < 1.3 from x 0.1-1% (e.g., |Martini et al. 
|2006| |2009l lEastman et aL| |2007l l. These luminous 
AGN fractions are approximately an order of magnitude 
lower than those found in field-galaxy regions over the 
same redshift. This decline in AGN activity with de- 
creasing redshift broadly tracks that found for the star- 
forming population in galaxy clusters (e.g., Saintonge 
leTaLl [20581 Finn et al.| |2ofol |Atlee et aL] |201 l| i, in- 
dicating a general shut-down of activity in the densest 
regions of the Universe over the last « 9 Gyrs of cosmic 
time. 

The results found for AGN activity in galaxy clus- 
ters suggest that growth of galaxies and BHs in over- 
dense regions must have been significantly more rapid 
at higher redshift. Indeed, we find direct support for 
this view from deep X-ray observations of the densest 
regions in the distant Universe (protoclusters). For ex- 
ample, in the z a; 3.1 SSA 22 protocluster, Lehmer et al. 
(2009b) found a factor of ^ 6 increase in the fraction 
of galaxies hosting AGN activity compared to the field- 



galaxy environment at z ^ 3.1, showing the opposite 
trend to that found for z < 1 .3 galaxy clusters (see also 
|Digby-North et aLl|2010[ ). Furthermore, there is tenta- 
tive evidence that the largest incidence of AGN activ- 
ity is in the densest regions of the protocluster dLehmer 



et al. 2009a I, suggesting that even the local environ- 
ment has a large effect on the triggering of AGN activ- 
ity, in contrast to AGNs identified in z < 1 overdense re- 
gions (e.g., Gilmour et al. 2007 Kocevski et al. 2009| l. 
The SSA 22 protocluster is predicted to become a large 
galaxy cluster by z ~ 0, with a similar m ass to that 
of Coma (a 10'^ Mq, [Steidel et al.||l998l l. The stark 



difference in activity between this z ~ 3 protocluster 
and z < 1.3 galaxy clusters of similar mass therefore 
suggests a dramatic change in fuelling mode, from a 
large abundance of cold gas available at high redshift 
to the predominantly hot gas when the galaxy cluster 
virialises. 

These studies show that environment has a major role 
in the growth of BHs. The optical and X-ray AGN activ- 
ity is likely driven by the availability of a cold-gas sup- 
ply. At high redshift this is most readily available in the 
densest environments but by lower redshifts the major- 
ity of the cold-gas in dense regions has been heated and 
cannot be accreted efficiently - at low redshifts AGN 
activity is more prevalent in typical regions of the Uni- 
verse. This shut down of activity is most dramatically 
seen in galaxy clusters and it may be related to the dark- 
matter halo. For example, the lack of significant AGN 
activity with dark-matter haloes » lO'^ Mq may be 
due to the halo virialising and shock heating the gas 



and suppressing vigorous mass accretion (e.g., Catta- 
|neo et al.| [2006,, 2008, Dekel & Birnboim, 2006, Keref 



et al. 2009 1. Indeed, luminous AGNs and quasars ap- 
pear to be universally hosted by dark-matter of 10'^- 
lO'^ Mq; see Fig. 5. The dark-matter halo may therefore 
control the cold-gas accretion and may strongly influ- 
ence the overall evolution of AGN activity with redshift. 

4. What fuels the rapid growth of the most massive 
(and also the first) black holes? 

In the previous section we discussed the host galax- 
ies and evolution of AGNs, which provide the overall 
picture for the joint cosmological growth of galaxies 
and their BHs. However, the work from the majority 
of these studies focused on typical AGNs with lumi- 
nosities near or below the break in the AGN luminosity 
function. While these objects represent the majority of 
AGN in the Universe and may play an important role in 
galaxy growth, they are unlikely to be responsible for 
the bulk of the accretion onto the most massive BHs. 



21 



Synthesis models of the BH population show that for 
systems with Mbh > 10^ Mq, the major growth phase 
occurs in very powerful, high-Eddington ratio quasars at 
high redshift (see §3.3.1). The large-scale process that 
is likely to triggers this rapid BH growth phase is galaxy 
major mergers, which is also expected to initiate lumi- 
nous star-formation activity. In this section we further 
discuss recent progress on the demographics of power- 
ful quasars and the physical mechanisms that fuel their 
rapid growth. We also explore the formation and evolu- 
tion of the initial BHs at high redshift that comprise the 
"seeds" which eventually grow into the BHs that power 
the quasars. 

4.1. Masses and evolution of rapidly growing black 
holes 

The foundations in understanding the rapid growth of 
massive BHs come from measurements of demograph- 
ics (BH masses, Eddington ratios, and space densities). 
In recent years the most powerful observation tool has 
been been spectroscopic studies of optically-selected 
quasars. Their high luminosity and characteristic col- 
ors mean that quasar luminosity functions can be traced 



out to high redshifts ( z < 7; fan et al.||2006[|Mortlock 
|et al.||20l"T] ), with the peak in the AGN space density of 
z * 1-3 (e.g., |Richards et al.|[2005||Croom et al!| [2009) 1; 
see §3.3.1. Compared to X-ray or IR observations, stud- 
ies of optical quasars are strongly biased against se- 
lecting obscured sources, but they have the key advan- 
tage of enabling detailed studies of the mass accretion 
(e.g., Floyd et al. 2009 Down et al. 2010| l and accurate 
BH masses (using virial mass estimators; e.g., Vester- 



gaard|2002 1 [Peterson et al.|2004[ [Vestergaard & Peter- 
son|2006t|Kollmeier et al.|20061|Shen et al.|20081|Me~ 
loni et al.|2010l l. 

The virial BH mass estimators are based on the sim- 
ple principle that the broad-line region (BLR) is in virial 
equiUbrium, such that 



BH - 



fv^RE 



(8) 



where / is a dimensionless quantity, of order unity, 
that depends on the geometry and kinematics of the 
BLRj^he velocity v can be measured from a single 
spectrum from the width of the broad lines, but the ra- 
dius Rblr is more challenging to estimate. All virial 
estimates are based on results from reverberation map- 
ping of local AGN, in which the time lag between the 
light curves of the continuum luminosity and the broad 
lines allows an estimate of Rblr (e g-, Kaspi et al. 2005 



[Bentz et"aLl |2009 ). These studies find a correlation of 
roughly R^lr <^ L'IL implying that Mbh v^L" ^ 



with a normalization that needs to be calibrated against 
other BH mass estimators and which depends on the 
kinematics and structure of the BLR. 

Repeated measurements of individual systems sug- 
gests that variation in the time lag and line widths re- 
sult in constant Mbh estimates, supporting the valid- 



ity of this technique (e.g., Peterson et al. 2004 Bentz 



[et al.| [2009| l. The local relation is generally calibrated 
on the width of Ha or HyS; however, at higher redshift we 
must use other lines (Mg n and C iv), which have also 



been validated from local measurements (e.g., McLure 



terson 



& Jarvis[ [2002{ [Vestergaard[ [2002) [Vestergaard & Pe 



2006 1 but for which the absolute caUbrations are 
quite uncertain (as discussed below). 

Despite these hmitations, a number of authors have 
forged ahead to produce Mbh estimates for many thou- 
sands of quasars from large spectroscopic surveys, most 



prominently SDSS (e.g., Shen et al.[ [2008^ 1 but also with 
AGES ([Koffmeier et al.| [2006|l, zCOSMOS ( |Merloni[ 



et al. 201 op and other surveys. Typical quasars have 
relatively large BHs (Mbh > 10^ Mq), with an efi'ec- 
tive upper limit around Mbh ~ 10'" Mq. Accurate BH 
mass estimates allow for a reasonable measurement of 
the Eddington ratio (see Eqn. 6 and Fig. 1). However, 
using optical spectra alone, these calculations require 
some estimate of the bolometric correction (BC) from 
the UV or optical continuum to obtain Lboi, which can 
depend on physical parameters (as discussed in § 2.3) 
such as luminosity, Eddington ratio, or BH mass (e.g.. 



Marconi et al. 




2004, Hopkins et al. 2007[ [Vasudevan 


& Fabian 2007 


Kelly etal. 2008 Davis &Laor 2011 


Raimundo et al. 


201 1 1. The Eddington ratio therefore 



becomes 



/I ccBC 



rO.5 



„2 ' 



(9) 



meaning that the observed distribution in A is primarily 
driven by the values of the broad emission line widths. 

These estimates generally produce a distribution in 
A which peaks around A = 0.1-0.3, declining rapidly 



2006 



for A > 1 and A < 0.01 (e.g., KoUmeier et al. 
Shen et al. 2008[ ); contrast with the Eddington ratios es 
timated for more typical X-ray AGNs (see §3.3.1). This 
confirms that optically-bright quasars do indeed repre- 
sent rapid growth phases, and so are in the regime of 



^For a spherical BLR geometry, / = 3/4 if the Hne width is ex- 
pressed in terms of the full-width half maximum (FWHM), as is com- 
monly adopted in virial BH mass estimates (e.g., Vestergaard 2002 
KoUmeier et al.| |2006l [Vestergaard & Peterson, 2006^ Shen et al. 



2008^ . Some studies have suggested, however, that FWHM is not a 
good way to characterize line widths, particularly for noisy data (e.g., 
[Fromerth & Melia|[2000| [Peterson et"ar][2004||Peterson[[20Tl) . 
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optically-thick, geometrically-thin disk accretion as dis- 
cussed in §2.3. The fact that few quasars are detected at 
A < 0.01 may result from the BH accretion flow switch- 
ing to the optic ally- thin mode at low Eddington ratios, 
but may also be the result of selection effects that would 
prevent such a low-Eddington source from having the 
blue continuum and broad lines characteristic of quasars 
Popkins etaL) [20091 ). 

Comparison of measurements at z ~ 6-7 to z < 
3 suggest that higher-redshift quasars are typically at 
lower Mbh but accreting at higher A than their lower- 



shown that a significant fraction (likely a majority) of 



redshift counterparts (e.g., Trakhtenbrot et al. 201 1 De 



|Rosa et ak] [201 l| l. However, interpretation of any ob 
served Eddington ratios is complicated by selection ef- 
fects that can strongly skew the A distribution. Indeed, 
a full Bayesian treatment of the observed distributions 
suggests that the typical intrinsic Eddington ratio for 
massive BHs at high z is closer to A ^ 0.05 or lower 
( |Kelly etal.||20T0l i. 

When interpreting these results it is crucial to de- 
termine the random and systematic uncertainties in the 
estimates of Mbh- Comparing Mbh measurements 
for nearby objects using multiple techniques (includ- 
ing stellar and gas dynamics from spectroscopy, anal- 
ysis of megamasers, and reverberation mapping), the 
error in virial Mbh estimators is found to be »0.4 dex 
(e.g. 



Onken et al. 



2004, Woo et al. 



est measurements of the radius-luminosity relationship 
suggest that its intrinsic scatter is only ^0. 1 dex (P eter-| 
[sori) 2010) . Thus in principle, radius-luminosity cali- 
brators should be highly effective, but there are impor- 
tant systematic uncertainties. One of these is the overall 
normalization of the virial estimator, which is related to 
the structure and kinematics of the broad-line emitting 
gas and is calibrated by comparison to estimates from 
the Mbh - cr relation (e.g., Onken et al. 



2004 Woo 



et al. 2010[ l. This cahbration is performed using nearby 



low or moderate-luminosity AGN, but at higher Edding- 
ton ratios radiation pressure can significantly change 
the kinematics and geometry of the BLR (e.g., |Mar-| 
Iconi et al.| [20081 INetzer & Marziani] [20T0l ). Further 
issues are that low signal-to-noise measurements may 



introduce large asymmetric uncertainties in Mbh (Den 



ney et al. 2009), and the value of luminosity included 



in the virial estimator may include a color term account- 



ing for the slope of the continuum ( Assef et al. 201 1 1. 
These effects may account for claimed discrepancies in 
the Mbh estimates using the UV and Balmer lines. 

Optically-bright spectroscopic quasars provide a 
wealth of information on demographics, but the sam- 
ples are strongly biased against AGN that are obscured 
by dust and gas. Recent multiwavelength surveys have 



luminous AGN at z > 1 are obscured (e.g., Tozzi et al. 



IHickox et all 120071 [Alexander et al.| |2008b 



Po zzi et al.| l20T0j [Bauer et al.| l20T0l [Treister et al 
20091 |2010bl. A complete census requires a com- 
bination of selection techniques, including searching 
for X-ray sources with extremely low ratios of opti- 
cal to X-ray flux (e.g., |Pozzi et all |2010[ |Bauer et al 



2010)), X-ray spectroscopy of known quasars (e.g.. Page 



etaTj |2004l [Comastri et~aLl 120111 ), detection of lu 
minous, narrow emission lines (e.g., 'Zakamska et al. 



20031 [Alexander et al.l l2008bl iReyes et al., ,2008, ,Vig 



nah et al.[ [20101 [Gifli et al.||2010| [Juneau et al.[[20TT]r 



or through optical-to-mid-IR SEDs (e.g., Hickox et al. 



2007 Bauer et al. 2010) . Since most obscured quasars 
do not show bright optical continua or broad emission 
lines, estimates of the bolometric luminosity and Mbh 
are more difficult. However, given these limitations, 
in general it appears that the luminosity and Edding- 
ton ratio distributions are broadly similar to unobscured 
quasars detected using similar techniques (e.g.,iHickox 
[eFaLl l2007j |Brusa et al.[ l20T0l [Mainieri et al.| |201 ij . 
Even with the most powerful contemporary techniques 
that employ a combination of deep IR, optical, and X- 
ray data, the most heavily obscured quasars (Compton- 
thick sources with A^h ^ lO^'* cm~^) remain extremely 



2010) . The lat- difficuh to identify (e.g.,lAlexander et"aLll2008b, Treis^ 



[ter et al.[[2010bt[Alexander et al.[[20rT] i7Therefore, esti- 
mates of the total obscured quasar population are either 
lower limits or depend on assumptions of the distribu- 
tion of absorbing column densities (A^h); however, these 
studies robustly demonstrate that a large fraction of the 
rapid BH growth is obscured. The nature of this ob- 
scured population may shed light on the processes that 
fuel the growth of the most massive BHs, as discussed 
below. 

4.2. Quasars, starbursts, mergers, and the evolutionary 
sequence 

Given a cosmological census of rapidly growing BHs, 
the next challenge is to understand the physical pro- 
cesses responsible for driving the accretion. The ac- 
cretion rates of up-to ^ 10-100 M© yr"' in quasars 
cannot easily be produced by secular processes that 
could trigger lower-luminosity AGNs (as discussed in 
§3.2.3 & §3.3.2). Powerful quasars require higher gas 
inflow rates and are more likely to be driven by gas- 



rich galaxy mergers (e.g., Kauffmann & Haehnelt 2000 



'Springel et al] [2005bt [Hopkins et al.','2006b'); however, 
large-scale secular instabilities may also be effective in 



particularly gas-rich distant galaxies (e.g.. Mo et al. 



23 




gas-rich 
galaxy(s) 



SMG/ULIRG 



obscured 
quasar 




unobscured 
quasar 



early-type 
galaxy 



Figure 6: Schematic diagram to illustrate the main components in the major-merger evolutionary scenario first proposed by |Sanders et al.|{l988^ . 



[T9981 [BowereTaLl [20061 IGenzel et al.| [20081 [Bourn 
[naud et aLl [201 Support for major-merger driven 
quasar activity comes from ( 1 ) the large fraction of sys- 
tems with disturbed morphologies (see §3.3.2), and (2) 
the good agreement between predictions for the merger 
rates from dark matter simulations (for adopted empiri- 
cal prescriptions the quasar fueling) and the spatial clus- 
tering and space densities of distant quasai-s (e.g., [Hop-| 
[kins et al.j [20081 [Treister et al.|[20T0ai ). 

While mergers are favoured by a number of mod- 



els (e.g.,[KaufFmann & Haehnelt 2000[ Springel et al. 


2005b 


Hopkins etaL 


|2006bi ISijacki et al. 2007 Di 


Matteo et al. 


2008|l, any quasar triggering mechanism 



requires relatively massive systems with large supplies 
of cold gas, which are generally f ound in dark mat ter ha- 
los with Mhaio ~ IO'2-IO'^ Mq ( [Croton[ Figure 4 [2009 
see also Fig. 5), just below the "maximal quenching' 
mass scales. Spatial clustering and environment mea- 
surements of quasars (e.g 



[eFaLl [20091 [Hickox et al.j [201 1[ [Carrera et al.| [201 Ij l" 



Ross et al. 20091 Lietzen 



suggest that quasars do indeed reside in halos of these 
masses at every redshift (as discussed in §3.4 and shown 
in Fig. 5). This implies that at high redshift, quasars are 
found in the largest collapsed system in the Universe 
(and are the progenitors of today's most massive early- 
type galaxies) while in the local Universe quasars are 
found in much more typical galaxy environments. Thus 
the mass of the dark matter halo may itself be the key 
parameter in understanding the fuelling of quasars. 

The rapid flow of cold gas that is necessary to fuel a 
quasar will inevitably be expected to also result in high 
rates of star formation (as discussed in § 2.1-2.2). Ro- 
bust evidence for links between powerful starbursts and 
quasars come from studies of local powerful IR galaxies 
(LiR > 5x10" Lq). The vast majority of such objects in 
the local Universe are major mergers of galaxies, with 
higher luminosities found during late stages when the 



galaxies are at small separations (e.g., Clements et al. 
[T996';'Ishidal. |2004| l. At higher Ljr the fraction of the lu- 
minosity from the AGN increases, and the large masses 
of nuclear gas and dust ensure that much of the BH 



growth is observed to be heavily obscured (e.g. , |Tran[ 
[eTaLl [200T1 [Yuan et all [MTOj [Iwasawa et al.j [20TTl 



Petric et al. 201 lb i. The local results are broadly con- 
sistent with models in which mergers fuel a rapid star- 
burst and a phase of obscured BH growth, followed by 
an unobscured phase after the gas is consumed or ex- 
pelled from the galaxy by stellar or quasar feedback 



(e.g., [Sanders et al.[[T988l[Di Matteo et all [20051 [H^ 



kins et al. 2008| l; see Fig. 6. However, powerful star- 
bursts are rare in the local Universe, compared to higher 
redshift where they dominate the star formation density 
(e.g., |Le Floc'hetaLl[2005l[Rodighiero etal.[|2010b| ). A 
key question, then, is whether a similar starburst-quasar 
scenario is the dominant process at high redshift, during 
the peak epoch of quasar activity where the largest BHs 
accreted most of their mass. Testing this picture is the 
subject of a number of recent studies. 

One approach is to select high-redshift starburst 
galaxies based on their IR or submmilimeter emission, 
and study the growth of BHs in these systems. Use- 
ful observational tools are X-ray observations and mid- 
IR spectroscopy, which can distinguish between dust 
heated by star formation and the AGN. The most pow- 
erful starbursts at high redshift, submillimeter galaxies 
(SMGs), have gas kinematics and morphologies that are 



characteristic of mergers (e.g., Tacconi et al. 2008 En- 
gel et al.[ [20T0l [Riechers et al.[ [201 l| l. A high frac- 



tion of these objects also host AGN (e.g., Alexander 
[eTall [2003a| [20051 [Laird et al.| fMO\ , but fliey gen- 
erally have Spitzer IRS mid-IR spectra that are domi- 
nated by star formation as indicated by luminous PAH 



emission features (e.g., Valiante et al. 2007 Pope et al. 



20081 [Coppin et all |20T0| . Only 15% of SMGs are 
dominated in the mid-IR by steep AGN continua, and 
even these powerful AGN generally do not produce the 
bulk of the bolometric output, which is dominated at 
longer wavelengths by the cool dust from star forma- 
tion (as also found for star-forming galaxies detected at 
70 jum; 'Sym eonidis et al.|20 10\ However, the presence 
of powerful AGN in some starbursts is consistent with 
a "transition" phase between powerful star formation 
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and rapid BH growth (e.g., Pope et al. 2008[ Coppin 
et al.||2010| as predicted by major merger fueling mod 



els (e.g.,|Springel et al.|[2005bl|Di Matteo et al.|[2005| 



Hopkins et al.||2008| l. In SMGs that host AGN activity, 
estimates of the galaxy luminosities and BH masses (as 
discussed in § 5.2; e.g., Alexander et al. 2008a[ Car- 



|rera et al.] |2011| ) indicate that the BH is undermassive 
relative to the host galaxy and may therefore be in an 
early stage of accretion that may precede the bright un- 
obscured phase. These results give broad support to a 
picture in which a rapid starburst is associated with the 
triggering of accretion onto a BH which quickly grows 
in mass. 

Further exploration of this "transition" population 
comes from studies of obscured quasars selected in- 
dependently of any star-formation signatures (such as 
far-IR or submillimeter emission). In the unified AGN 
model, obscuration in quasars is purely due to orienta- 
tion and so the host galaxies and large-scale structures 
of quasars would be independent of obscuration (e.g., 
Antonucci] [T993f In contrast, evolutionary scenarios 



would suggest significant diff'erences in the host-galaxy 
proper ties (e.g., [S anders et al. 1988; Springel et al.j 
|2005b|[DrM^eo et al., 2005,. Hopkins et al. , 2008) ; see 
Fig. 6. Some studies have focused on the SEDs of X-ray 
selected obscured quasars, which can often be modeled 
with an extension of the simple "torus" geometry, but 
in many cases also require a contribution from the host 
galaxy at longer wavelengths indicating an energetically 
important st arburst (e.g., Mainieri et al., 2005a; Vignali| 



^FaLl [20091 iPoUetta et al.[|2008[|Pozzi et al.[|2010| l. 

Additional tests of this evolutionary scenario come 
from clustering analyses, as diff'erences in host halo 
masses would be inconsistent with a simple "unified" 
picture in which obscuration is purely due to the orien- 
tation of an obscuring torus along the line of sight (e.g., 
Antonuccij [T993| |Urry & Padovaiiil [1995) . The first 
measurement of the spatial correlation of IR-selected 



obscured quasars (Hickox et al. 20111 indicates that 
they are clustered at least as strongly as their unobscured 
counterparts, with Mhaio possibly above lO'"^ /i"' Mg. 
This suggests the possibility that obscured quasars rep- 
resent an earlier evolutionary phase in which the BH 
is undermassive relative to its host halo, such that ob- 
scured quasars are found in more massive halos com- 
pared to unobscured quasars with similar BH mas^ 

'We note that lower-luminosity X-ray selected AGN appear to 
show the opposite trend, with stronger clustering for unobscured com- 
pared to obscured sources jAIlevato et al.||2"01H . This suggests that 
this evolutionary scenario may not hold for lower-luminosity sys- 
tems, and provides strong motivation for more precise future cluster- 
ing measurements that probe a wide range of luminosity and redshift. 



If obscured quasars represent an early phase of ac- 
cretion then do we see the shutdown of star formation 
upon the emergence of an unobscured quasar? Evidence 
based on the submillimetre properties of quasars sug- 
gests that unobscured quasars do coincide with a signifi- 
cant decrease in star formation. |Page et al. (2001 2004[ l 



and [Stevens et al.| ( |2005| l find a stark difl'erence in the 
submillimetre emission between X-ray obscured optical 
quasars and X-ray unobscured optical quasars: the X- 
ray obscured optical quasars (which represent x 15% of 
the optical quasar population) have SFRs implied from 
the submillimetre data that are about an order of magni- 
tude higher than the X-ray unobscured optical quasars. 
Tentative evidence for the catalyst of this decreased 
star formation is found from the spectra of these sys- 
tems. The X-ray spectra of the X-ray obscured optical 
quasars indicate that the obscuration is due to an ionized 
wind, providing evidence for a quasar-driven outflow 
that may be associated with the short-lived transition 



from starburst to AGN-dominated systems (Page et al. 



201 1) . Near-IR IFU spectroscopy of another X-ray ob- 



scured optical quasar shows evidence for an energetic 
outflow on scales of » 4-8 kpc, which may be starting 
to shut down the star formation across the whole galaxy 
( [Alexander et al.| [2010| l. At lower redshifts, massive 
and energetic molecular outflows have been observed in 
powerful AGN (Sturm et al. 201 l| l including the broad- 



Tllli 

absorption-line quasar Mrk 231 ( [Fischer et al. 2010|l 



for which the outflow is extended on kpc scales ( Fer-[ 



uglio et al. 2010 1. These results provide tentative evi- 



dence that AGN-driven winds may indeed influence gas 
on the scale of the host galaxy and help quench the for- 
mation of stars (as discussed in § 5). 

Multiwavelength studies of optical quasars also in- 
dicate a small population 10%) of objects that lack 
clear hot-dust signatures (e.g., Jiang et al. 2010 Hao 
et al.| [2010 201 1\ . There are a number of interpre- 
tations for the lack of hot-dust emission in these ob- 
jectsj^ncluding the possibility that some hot dust-poor 
quasars represent the end stage of the quasar evolution- 
ary sequence. 

Taken together, the populations of the most-active 
systems at z x 1-3 may be tentatively placed 
into a broad evolutionary sequence (SMGs/ULIRGs- 
obscured quasars-unobscured quasars-"hot dust-poor" 
systems); see Fig. 6. However the precise relative num- 
ber densities, duty cycles, and environments of these 
diff'erent objects remain poorly constrained, and it is 



'"in the highest-redshift (z > 6) systems, there may simply not be 
enough dust in the Universe to form a significant dusty torus jjiang| 



[etaT] 



nmo). 
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therefore not clear how well the different populations 
can be accounted for by a single evolutionary scenario. 
Future observations with larger samples and improved 
diagnostics will provide the tools to verify or falsify this 
general picture for the rapid growth of distant massive 
BHs. 

4.3. Very high redshifts and the formation of "seed" 
black holes 

Through the preceeding sections we have discussed 
the significant progress made in understanding the pro- 
cesses that trigger the rapid growth of the most mas- 
sive BHs. However, any scenario of quasar fuelling re- 
quires an initial BH onto which accretion-driven growth 
can occur, and the nature of how and when these initial 
seeds were formed remains an important open question. 
In many cosmological simulations that include BHs, 
these "seed" BHs are simply put in by hand with some 



arbitrary mass (e.g.. 


Bower et al. 2006 2008 Booth & 


Schaye 2010 Panic 


lakisetal. 2011 2012|l. The late- 



time properties of the BHs and galaxies are generally 
insensitive to the mass of the seeds, as later accretion 
"erases" any memory of the initial conditions. However, 
the precise nature of the seed BHs is still crucially im- 
portant, as the early formation and growth of BHs may 
in fact represent a key component in the early formation 
of structure. Early BH accretion may also produce a 
significant fraction of the total radiation background re- 
sponsible for reionization of the Universe at z i 7 (e.g., 
IMadau et al.| [2004, MiraberetaL] [2011] ). 

One primary challenge for any model of seed BH for- 
mation is posed by the existence of luminous quasars 
with Mbh > 10^ Mq at z > 6, when the Universe was 



<1 Gyr old (e.g.. Fan et al. 2006 Jiang et al. 2009 



jWillott et al.1 '2010 ; De Rosa et al. , 201 1 ; Mortlock etaL 
[201 \\. For a BH accreting with an Eddington ratio A and 



assuming a radiative efficiency e, the growth time of a 
BH from an initial mass M,„ to a final mass My,„ is given 



by (e.g., [Volonteri,.2010j : 



fgrowth = 0.45 Gyr A ' In 

1 - e 



M 



fin 



(10) 



This formation time scale puts robust lower limits on 
the masses of the seed BHs, for there must be sufficient 
time in the history of the Universe for them to grow 
into the massive BHs seen in the high-z quasars. For 
an initial BH formed at a high redshift (say z - 20) 
to reach 10^ a;0.7 Gyr later at z = 6, assuming 
standard e = 0.1 and continuous accretion at the Ed- 
dington limit, it must have started with M,>, ~ 100 Mq 
(with the caveat that some additional growth can oc- 
cur through BH-BH mergers; e.g., Sesana et al.[[2007" 



[Arun et al.[[2009t[Sereno et al.|2011| l. For a lower for- 
mation redshift or more stochastic accretion (and thus 
lower average A) intial BH masses would necessarily 
have been higher. A successful model of seed forma- 
tion must therefore produce sufficiently massive BHs to 
satisfy these constraints. 

There are currently three main candidate mechanisms 
for seed formation: (1) remnants of massive population 
III stars (2) direct collapse of primordial gas clouds, and 
(3) runaway collisions in dense stellar clusters. Each 
model makes different predictions for the masses and 
number densities of the BH seeds (see [ Volonterij [20 lOl 
for a detailed review). 

The first candidate mechanism proposes that BH 
seeds are produced by the deaths of Population III 
stars. Stellar models suggest that this first generation 
of stars, which formed from very low-metallicity gas, 
can have large masses of ^ 100-600 Mq (Abel et al. 
[2002[[Bromm et al.|[2009| ), and that the collapse of stars 
with masses > 260 Mq can produce ~100 Mg BH rem- 



nants (e.g.. Fryer & Kalogera 2001 Madau & Rees 



2001[ l. This model has the advantage of following the 



well-established mechanism for the formation of lower- 
mass BHs through stellar processes, but is limited by 
the small mass of the seeds. Radiative feedback limits 
the possible accretion rates onto seed stellar-mass BHs 
(e.g., Johnson & Bromm 2007 Alvarez et al. 2009| l, 



making it difficult to construct mechanisms by which 
such BH seeds can achive the high accretion rates re- 
quired to produce the high-mass quasars observed at 
z > 6 (e.g., Haiman & Loeb 2001 Volonteri & Rees 



2005 Volonteri et al. 2006| l. Therefore, stellar pro 
cesses may be responsible for a significant population 
of early BHs, but in the most biased regions (where we 
find the high-z quasars) another formation mechanism is 
required. However, recent studies have challenged the 
idea that the first stars are very massive, due to either 
fragmentation or feedback effects (e.g., 'Glover et al 



,2008 , McKee & Tan , 20081 [Turk et al., .2009, .'^enti^ 
[Stiavellil [2009) [Stacy et al.| [20T0l [Greif et al.| [20n] i, 
casting some doubt on the the role of Population III rem- 
nants as BH seeds. 

The second candidate mechanism proposes that BHs 
form directly inside the densest peaks in the matter dis- 
tribution. In certain circumstances, the BHs can be pro- 
duced through the coUapse of gas clouds. In regions 
with radiation fields that can dissociate Ho molecules, 
star formation is suppressed and gas cooling proceeds 
through atomic processes. Detailed hydrodynamic sim- 



ulations have explored this process (e.g., |Wise et al 
20081 [Regan & Haehnelt| [20091 [Shang et al., .2010 



Johnson et al. 201 1 1 and indicate that supermassive 
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stars at the centers of halos can, in principle, accrete up 
to a few xlO'' Mq over ~2 Myr (John son et al. 2011 1. A 
significant fraction of this mass may then go into a seed 
BH which can continue to grow, although less rapidly 
because of the effects of radiative feedback. A related 
possibility is that direct collapse does not occur in an 
isolated halo, but is triggered by the mergers of multiple 
halos, providing highly unstable c onditions and large 
gas inflow rates > 1 Mq yr ' (e.g., Volonteri & Begel- 
2010[ Begelman 2010 1. Such conditions can pro 



duce "quasistars", massive optically-thick structures in 
which the radiation pressure from accretion is balanced 
by the ram pressure of infalling gas, in which gas accre- 
tion can significantly exceed the Eddington limit for the 
central BH ( [Begelman et al. 2006[ l. This process can 
produce BH seeds with masses of up to ~10^ Mq in the 
most biased regions of the early Universe. 

The third candidate mechanism proposes that the 
seed BH formation is preceded by the formation of 
dense stellar clusters, the centers of which then un- 
dergo runaway merging due to stellar-dynamical pro- 



cesses that produces a BH (e.g., Devecchi & Volonteri 



[20091 Pevecchi et al.| |20T0l [Davies et al.| |201 l| l. This 
process can proceed in regions with significant molecu- 
lar gas, and does not require extremely high inflow rates 
as in direct gas-dynamical collapse, but generally pro- 
duces smaller BH seeds with masses ~10^ Mq. 

Given these possible mechanisms for the formation 
and growth of seed BHs, the challenge for observers 
is to distinguish between them. One possibility is to 
trace the high-z BH population through accurate mea- 
surements of the number density of AGN at z > 4. The 
most progress has been made through studies of high-z 



optical quasars (e.g.. Fan et al. 2006 Jiang et al. 2009 



[Willott et al.l|20T0 ';'De Ros a et al.i|201 1 [[Mortlock et" 
[201 l |l. However, these samples comprise only the most 
massive end of the BH distribution, and these quasars 
appear to be growing more rapidly than their lower- 
redshift counterparts ( jPe Rosa et al. 2011 1 indicating 
that they are also biased to higher Eddington ratios. In 
principle. X-ray observations allow us to explore further 
down the mass distribution, and some success has been 
achieved by searching for X-ray counterparts to high- 
redshift optical and IR sources (Fiore 2010[ ). In the 



future, a combination of sensitive, high angular resolu- 



tion X-ray surveys (e.g., eR OSITA,\Predehl et al [2007| 
Wide-Field X-ray Telescope, 'Murra y et al.|2010 [Brusa 



et al.[[20lT I over wide areas covered by IR and optical 



observations will be required to accurately measure the 
high-z AGN luminosity function to faint limits. 

Measurements of the AGN space density are sensi- 
tive not only to the seed formation mechanisms, but 



also the rate of BH growth and BH-BH mergers. In 
principle it may be possible to observe the seed BHs 
directly, by searching for the signatures of accretion 
onto < 10^ Mq BHs. Even for such BHs, the hard 
radiation field may produce characteristic nebular line 
emission (in particular. He II /1 1640) that could be de- 
tected with NIRSpec on the James Webb Space Tele- 
scope Pohnson et al.[[20TT] l, while quasistars could pro- 
duce a significant population of sources observable at 



mid-IR wavelengths with JWST ( Volonteri & Begelman 



2010i In the more distant future, these BHs may also be 



detectable through gravitational wave signatures when 
they merge (e.g.,|Sesana et al.[ [2007) [Arun et al.| [20091 



[Sereno et all [201 1| ). 

Finally, a complementary strategy for understanding 
seed BHs is to look for residual signatures at lower red- 
shifts. For the most massive systems, the properties of 
the seeds have been "washed out" by subsequent mass 
accretion events. However, some of the lower- mass 
(< 10^ Mq) BHs in the local Universe may represent 
relics of the initial epoch of seed formation (e.g., |van[ 



Wassenhove et al. 2010 1. Therefore the formation pro- 



cess for seed BHs may be reflected in the distribution of 
low-mass BHs within dark matter halos at low redshift, 
which can be probed using accurate measurements of 
the relation between BH masses and galaxy properties 
(e.g., [Volonteri et al.[ [2008) . By studying seed BHs 
in concert with the subsequent fueling and evolution of 
powerful quasars, we can aim to produce a coherent pic- 
ture of the growth of massive BHs from their early ori- 
gins to the present epoch. 



5. What is the detailed nature of AGN feedback and 
its effects on black-hole fueUing and star forma- 
tion? 

A recurring theme in any discussion of BH growth is 
the impact that energy released by accretion has on the 
surrounding gas in the system. In principle, a growing 
BH releases plenty of energy to impact its surroundings: 
the thermal energy of the hot gas atmosphere in a lO'"^ 
Mq dark-matter halo is «ilO^' erg, while the total accre- 
tion energy of a 10^ Mq BH is ^ 2x 10^^ erg! Therefore, 
given a sufficiently strong coupling between the radia- 
tive or mechanical output of the BH and the surround- 
ing gas, the AGN should be able to disrupt its environ- 
ment and potentially regulate its own growth and star 



formation in the host galaxy (e.g. Silk & Rees 1998 



Bower et al. 2006 1. In this section we explore the vari- 
ous physical processes, such as winds and jets, by which 
energy from the BH can couple to the surrounding gas. 
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We also address the observed consequences of this feed- 
back, with emphasis on the relations between BHs and 
their host stellar bulges, as well as the gas content of 
dark-matter halos. 

5.7. Driving feedback through winds and jets 

Any feedback process from a growing BH requires 
coupling between the energy released by the BH and 
the surrounding matter In general these can take two 
forms: (1) "winds" (often referred to as superwind- 
mode), which comprise wide-angle, sub-relativistic out- 
flows and tend to be driven by the radiative output of the 
AGN, and (2) "jets" (often referrred to as radio-mode), 
which are relativistic outflows with narrow opening an- 
gles that are launched directly from the accretion flow 
itself; see Fig. 7. As described in §3.32, §3.3.2, & §4.2 
the radiatively-dominated AGN that drive winds are ex- 
pected to be relatively high-Eddington ratio systems, 
while jets are most commonly produced (except for the 
highest-power sources) by lower-Eddington ratio accre- 
tion flows. 

5.1.1. Radiatively-driven winds 

In recent years, a wide range of observational ef- 
forts have been aimed at searching for evidence for 
AGN winds, largely through the presence of highly 
blueshifted absorption and emission lines. Much work 
has focused on bright, nearby AGN, for which high 
signal-to-noise observations readily allow the detection 
of the relevant line features. As discussed in §2.1, 
IFU observations of local Seyferts show evidence for 
outflowing gas on scales of ~10-100 pc (e.g.. 



"Superwind" 

expels cold gas reservoir 



"Radio" 

reheats cooling atmosphere 



Bergmann et al. 



IDavies et al. 2009 



Storchi- 



Storchi- 



Bergmann eTaL] |20To| [Schnorr Miiller et al.| |2011| ), 



while spectropolarimetry of a low-redshift quasar shows 
a high-velocity (~4,000 km s ') outflow within the nu- 
clear torus itself, close to the accretion disk (Y oung| 
2007| l. High- velocity winds are also commonly 



et al. 



observed using X-ray spectroscopy, which can trace 
atomic transitions in ionized gas. The presence of 
highly blueshifted absorption and emission lines in a 
number of quasars and Seyferts (e.g., PDS 456, PG 
I211-H143, NGC 4051, 3C 455, MR2251-178) indi- 
cate outflows of V a! 0.1-0. 3c ( [Pounds et al. | |2003 



[Reeves et aLl|2003i|Tbmbesi et al.|[20T0l|Lobban eTaL 



2011 Gofford et al. 2011 1, as expected for models of 
momentum-driven winds, as discussed below. A lower 
velocity outflow (~500 km s"') is observed in the local 
Seyfert 2 galaxy NGC 1068 (although it is confined to 
the nuclear regions rather than extending on host galaxy 
scales). The high-resolution X-ray spectrum for NGC 
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Figure 7: Schematic diagrams to illustrate the two main modes of 
AGN outflows: "superwind"-mode outflows such as those found in 
luminous AGNs and "radio"-mode outflows such as those found in 
low-excitation radio-loud AGNs. 



1068 is characteristic of photoionization rather than col- 
lisional ionization, suggesting that the outflow is driven 



by a radiatively efficient wind rather than a jet (Evans 
[eFaLllIOTOl i. 

Similar statistical analyses show evidence for power- 
ful outflows in distant quasars. In the spectra of high-z 
quasars from SDSS, the peak of the broad C iv emission 
is typically b lueshifted, indicating outflo ws with ~1000 
km s ' (e.g., Richards et al. 2002 201 1 1. Quasars with 



larger blueshifts show lower equivalent widths of C iv 
as well as weaker observed X-ray emission. A similar 
trend with X-ray emission is found for broad absorp- 
tion line (BAL) quasars, which are believed to represent 
the 20^0% of quasars (with evidence for a redshift 
dependence; e.g., Allen et ar|[2011 1 for which the nu- 
clear continuum is viewed through the outflow (there is 
evidence that maybe all quasars host energetic winds; 
Ganguly & Brotherton|2008] l. The maximum blueshift 



in BAL quasars (corresponding to the terminal veloc- 
ity of the wind) is higher for sources with weaker X- 



ray emission (e.g., Gallagher et al. 2006 Gibson et al 
[2009| l. Both these and the results on C iv blueshifts 
can be understood in terms of radiation driving of the 
quasar wind; with a weaker X-ray continuum (possibly 
associated with an absorber near the midplane), the gas 
close to the BH is less highly ionized and so more easily 
driven by the UV continuum radiation from the quasars. 

Evidence for outflows also comes from statistical 
studies of large numbers of moderate-luminosity AGNs 
at z i 0.4. A number of narrow-line Seyfert galaxies 
show complex [O in] emission with a broad blue wing in 
addition to the main narrow component (e.g., [Veilleux[ 
1991 Boroson 2005 | l, and statistical studies suggest 
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these may be ubiquitous (MuUaney et al., in prep). The 
width and luminosity of the broad component increase 
with L[oiii] but are independent of radio-loudness, possi- 
bly indicating that the winds are driven by the radiative 
output of the AGN rather than by relativistic jets. 

These studies demonstrate that energetic winds are 
common in the AGN population. However, they cannot 
establish on the basis of this data whether these winds 
generally have significant effects on the scale of the host 
galaxy since, in the majority of the cases the winds are 
observed only along the line of sight and there are no di- 
rect constraints on the spatial distribution of the outflow- 
ing gas (e.g., [Tremonti eFall |2007| |Arav et al.| [20081 
Dunn et al. 2010). To fully understand the mass distri- 
bution, velocity structure, and energetics of these winds 
requires improved constraints on the covering fraction 
and dumpiness the outflows, which are difficult or im- 
possible to determine from spectroscopy along a single 
line of sight. 

Recent observations have begun to overcome some 
of these challenges through spatially resolved spec- 
troscopy (IFU observations). For example, a quasar 
hosted by a z ~ 2 SMG shows broad blue wings to [O m] 
lines extend out to several kpc ( [Alexander et al.| [2010[ ). 
The wind velocity in this case is a; 300 km s"', which 
suggests that it could possibly be driven by the power- 
ful starburst instead of the AGN. In addition, the ob- 
served outflow is sufficient to strongly disturb gas in the 
galaxy, but is unlikely to completely unbind the mate- 
rial from the system (see §5.1.2 for IFU observations of 
several high-z radio-loud quasars). At lower redshifts, 
long-slit spectroscopy of luminous narrow-line quasars 
at z < 0.5 has revealed ubiquitous galaxy-scale ionized 
clouds with high velocity dispersions, which appear to 
be powered by the AGN and may be associated with 
outflows (Greene et al. 201 l| l. An extreme example of 
an outflow has also been observed with CO line obser- 
vations of the nearby quasar Mrk 231, in which broad 
wings to the observed lines indicate the presence of a 
high- velocity (~750 km s ') wind that is resolved on 
kpc scales and contains ~700 Mg yr ' of molecular gas 
( [Feruglio et al. 2010^). This outflow may be expected 
to evacuate the cold gas content of the galaxy within 
~10 Myr. Observations are revealing more and more 
examples of powerful feedback from radiatively driven 
winds, but the overall prevalence of these outflows and 
their overall impact on the galaxy population remains to 
be determined. 

5.7.2. Relativistic jets 

In contrast to the picture for radiative winds, relativis- 
tic jets from AGN are commonly observed to influence 



gas on the scales of the host galaxy or even its parent 
dark matter halo. Indeed, the brightest structures ob- 
served in radio AGN are often on > kpc scales and are 
produced by the coupling of the AGN outflow to its en- 
vironment. The physical parameters that produce AGN 
jets (low accretion rates, fueling from hot atmospheres, 
BH spin) remain poorly understood; however, the im- 
pact of the jets on the hot gas in galaxies and dark-matter 
halos is relatively well established. 

In massive systems with hot ionized atmospheres, 
AGN jets interact most strongly with the gas by in- 
flating bubbles of relativistic plasma in the ambient 
medium. In the classical picture for radio lobe forma- 



tion ( Scheuer 1974), the bubbles are strongly overpres- 
sured and so produce shocks that can directly heat the 
atmosphere. However, Chandra observations of hot gas 
around AGN jets yield only a few clear examples of 
large-scale shocks, and these are predominantly in low- 
luminosity AGN with FR I morphologies (e.g., [Forman[ 



et al. 2007 Croston et al. 2009| , as opposed to the 



powerful FR II sources that are expected to drive the 
strongest shocks. Therefore most of the mechanical en- 
ergy from the jets is transferred in the form of PdV work 
by the inflating bubbles (e.g.,[Birzan et al.[[2004[[Birzan[ 



etal. 2008| Cavagnolo et al. 2010 1 rather than flirough 



strong shocks. The bubbles are observed to carve out 
cavities in the hot gas, so there is expected to be little 
mixing between the hot gas and the relativistic plasma. 
In this case the pressure is simply the external pressure 
Pext (which can be measured from the X-ray images), 
so that AE - Pext^iobe- Assuming a timescale for the 
lifetime of the radio source (for example from the buoy- 
ant rise time of the inflated bubble), we can then es- 
timate the total mechanical power f mech output by the 
AGN (e.g., Birzan et al. '200T; Dunn et al. "SOPS'; Dunn [ 
CTabian, ,2008, ,Birzan et al., ,2008, ,Cavagnolo et alf 



2OTOI 1. 



Such estimates show that for most radio-loud AGN, 
the mechanical power in the relativistic outflows vastly 
exceeds (by as much as ~1000 times) the radiative out- 



put of the AGN (e.g., Birzan et al. 2008 Cavagnolo 



[et al.| [2010[ ). These studies also find correlations be- 
tween the radio luminosity and (much larger) mechani- 
cal power, which can be used to estimate overall distri- 
bution of mechanical energy input in the Universe (e.g., 
IHeinz et al. j [20071 [Merloni & Heinz| [2008i |La Franca[ 
et al. 2010 1; see §3.3.1. Analysis of the cavities can 



also be used to explore whether the mechanical output 
of the AGN is sufficient to balance the energy lost in 
the atmosphere through radiative cooling. While limited 
mainly to bright, massive clusters (where the cavities 
are most easily identified), the PdV work done by the 
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cavities can usually offset radiative losses (e.g., Brrzan 
etaL] [2004'; Dunn & Fabian 2008V Indeed, the heating 



from intermittent AGN activity and the radiative cool 
ing are often in relatively close balance ( [Rafferty et al. 



2008 1. These results suggest a limit cycle in which the 



cooling of hot gas onto the nucleus of the galaxy fuels 
intermittent AGN outbursts, which in turn provide suf- 
ficient heating to slow down or stop the accretion flow 



(e.g.. Best et al. 2005 Pope 2011) . This scenario rep- 
resents "feedback" in the true sense of the term, in that 
the energy output by the AGN has a direct impact on its 
own fuel supply. 

While most systems are roughly consistent with 
obeying such a limit cycle, there are some examples 
in which the radio jets impart far more power than is 
necessary to simply reheat the gas, but instead may be 
sufficient to unbind much of the hot atmosphere. Such 



sources are common among galaxy groups (Giodini 



|et al.|[20T0l ), but the most extreme case is in the massive 



cluster MS0735.6-H74 21 (e.g., [_M cNamar aet aL| [2005: 



[2009 j ; see also ,Siemiginowska et al. (2010 i. This ob- 
ject exhibits enormous cavities of ~200 kpc in diame- 
ter, for which the PdV work required is ~10^^ erg, or 
close to the entire accretion energy of a lO** Mq BH. 
The AGN in this source may be driven by the accre- 
tion of cold gas rather than fueling directly from the 
hot atmosphere, and may be an example of a radio-loud 
quasar (that is, a source with both powerful radiative 
and mechanical output), rather than lower-Eddington, 
optically faint AGN. In support of this picture, high- 
excitation (radiatively efficient) radio AGN are gener- 
ally found in blue star-forming galaxies that are ex- 
pected to be rich in cold gas, while low-excitation (ra- 
diatively inefficient) sources are found in red, passive 
galaxies with low cold gas content (e.g.,'Smolcic','2 00"9{ 
[Herbert et al] [SoTOj [Smolcic & Riechers, ,2011) . In 
the case of MS0735.6H-7421, the enormous mechanical 
power needed to produce the cavities has prompted a 
more extreme possibility, that the outflow must be pow- 
ered not only by accretion energy but also by the extrac- 



tion of BH spin energy by the relativistic jet (McNamara 



etal. 2009 2011 1. Whatever the ultimate power source, 



these extremely powerful mechanical outflows may be 
responsible for expelling significant gas from the cen- 
ters of halos, as discussed in § |5.2.3[ 

FinaUy, there are examples of relativistic jets interact- 
ing not only with hot, virialized plasma, but also with 
cooler gas in the host galaxy. IFU observations of pow- 
erful ( i 10^^ W Hz ' at 1.4 GHz) radio-loud quasars 
have revealed outflows in [O m] that are aligned with 
the radio jet axis over scales of ~10 kp c, with v ~ 1000 
andPWHM ~ 1000 km s ' (e.g.[Nesvadba etal. 



2006|[2007[[2008l ). These results suggest that relativis- 
tic jets can ionize and expel significant amounts of cool 
(^lO** K) gas, in addition to strong interactions with the 
virialized hot atmosphere. 

5.2. Observational consequences of feedback 

The discovery of the ubiquity of massive BHs in 
galactic centers and the apparently tight relationship be- 
tween the mass of a central BH and the properties of 
its host stellar bulge has motivated a great deal of inter- 
est in the physical consequences of BH growth. These 
results provide compelling (although indirect) evidence 
for interplay between the growth of the BH and the for- 
mation of the host galaxy, despite enormous differences 
in mass 3 orders of magnitude) and linear scale (^ 8- 
9 orders of magnitude). Given the energetic importance 
of BH accretion in the context of the galaxy, feedback 
from AGN was immediately recognized as a possible 
mechanism to produce these observed correlations (e.g.. 



Silk & Rees T998| Fabian 1999 King 2003 1 



5.2.1. Correlations between BHs and galaxy spheroids 
The first step in understanding BH-spheroid relation- 
ships is to measure accurately the nature and intrinsic 
scatter in the relations, as well as their evolution with 
redshift. A wide range of studies have explored correla- 
tions of Mbh with bulge mass, luminosity, velocity dis- 
persion, and gravitational potential, and found a remark- 
ably small intrinsic scatter for massive galaxies in the 
local Universe (e.g., Magorrian et al. 1998 [Gebhardt 



leFaLl [20001 IFerrarese & Merrittj |2000[ |Giiltekin eFal 
[2009| l. Recent studies of lower-mass galaxies suggest a 
change in relationship at the low end, with an increase in 
scatter (e.g., Greene et al. 2008 2010"), although some 
of this may be driven by the prevalence of pseudobulges 
in disk galaxies rather than classical bulges, where the 
correlations are much weaker or non existent (e.g., [Ko 
rmendy & Bender 201 1, Jiang et al. 2011 



To understand the physical origin of these corre- 
lations, important clues come from the evolution of 
the BH-spheroid mass relationship with cosmic time, 
which can now be explored observationally. The diffi- 
culty in measuring BH masses for distant galaxies re- 
stricts the current studies to systems hosting AGN ac- 
tivity. BH mass estimates are computed from either the 



km s 



virial method (for broad-line sources; e.g.. Woo et al. 
,20081 |Kim et al.||20081 [D ecarU et al.|[20i0b, Bennert 
[et al.|2010| [Merloni et al.|2010| see §4.1) or via some 
assumptions about the relationship between obscured 
and unobscured AGNs (as for SMGs; [Alexander et al. 
2008a| l, each of which have significant uncertainty. The 
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presence of an AGN also presents the additional chal- 
lenge of separating the light from the host galaxy to that 
from the accreting BH, and difficulties in separating the 
mass or potential of the stellar bulge from the rest of the 
galaxy (e.g., |Decarh et al.||2010"al [Merloni et al.||2010[ ). 

Despite these challenges, a number of authors have 
explored the evolution in MsH-Mbuige for distant AGNs. 
Studies of broad-line quasars consistently find that the 
BH is overmassive relative to its host galaxy compared 



to the local relation (e.g., [Walter et al. 2004 


|Woo et al. 


12008 ; Decarli et al.||2010b| Bennert et al. 


20101 Mer-I 


|loni et al.,_2010i. In most cases comparisons are made 



between Mbh and the total galaxy mass, rather than the 
spheroid properties, which are challenging to determine 
with current observations. These results have prompted 
speculation that BH growth precedes that of the host 
galaxy. However the observed evolution may be due in 
part to selection effects that cause more massive BHs to 
be more readily detected (e.g., [Lauer et aL| |2007[ |Shen| 



& KeUy 2010 1, and some observational results have 
shown httle redshift evolution, at least in the relation 
between BH and total galaxy mass (e.g., jJahnke et al. 



2009| l. In contrast, corresponding analyses of X-ray se- 



lected AGN in powerful starbursts (SMGs) indicate that 
the BHs lie up-to an order of magnitude below the re- 
lationship for z ~2 quasars, suggesting that these BHs 
are in the processes of "catching up" to their final mass 
( [Alexander et al.||2008a||Carrera et aLl|201 l| l; however, 
SMGs may represent an earlier stage in the evolution of 
quasars (see §4.2 and Fig. 6). At present, the signifi- 
cant uncertainties in the current approaches, as well as 
the biases inherent in limiting studies to galaxies host- 
ing AGN, make it difficult to draw robust conclusions 
about the evolution of BH-spheroid mass relationships. 
Improved observations and larger, well-characterized 
galaxy and AGN samples should provide better con- 
straints on this evolution in the future. 

5.2.2. Physical drivers of BH-spheroid relationships 

Redshift evolution notwithstanding, a great deal of 
theoretical eff'ort has attempted to explain the observed 
BH-spheroid relationships in the local Universe, with a 
particular focus on the effects of outflows from AGN. 
As one example of the vast array of analytical models. 



King et al. (2011 1 consider roughly Eddington-limited 
accretion in which the optical depth to photon scattering 
is ~1. In this case, the outflow momentum is compa- 
rable to the photon momentum, which yields a typical 
outflow velocity v ^ el Ac ~ 0.1c, as observed in the 
X-ray spectra of many AGN (discussed below). Comp- 
ton cooling by the AGN radiation field allows this out- 
flow to cool rapidly inside a radius of ~ 1 kpc, reducing 



the thermal energy of the flow but conserving its mo- 
mentum. Setting the thrust from such a "momentum- 
driven" flow to balance the gravitational force on the 
bulge gas naturally produces a BH-spher oid mass re - 
lationship with Mbh cr'^, as observed (King 2005|l. 
This model also predicts that BHs in AGN should in 
general be accreting near the Eddington limit, and grow- 
ing up toward the Mbh-o" relation that is observed for 



passive systems ( King 2010a) . 

Such analytical models provide a clear physical pic- 
ture for simple generalized conditions (spherical sym- 
metry, smooth matter distributions, and constant ac- 
cretion rates). An alternative approach is to study the 
growth of BHs in the context of large cosmological 
hydrodynamic simulations, which provide more real- 
istic physical conditions but have the hmitation that 
they do not directly resolve the sphere of influence of 
the BH. Instead, a common approach is to assume that 
some fraction of the gas inflowing into the central ~kpc 
(through mergers or other processes) accretes onto the 
BH, up to the Eddington limit (e.g., 'Springel et al. 



2005b Pi Matteo et al. 2005 2008; Booth & Schaye 



2010 1. This prescription is broadly valid if the transport 
of gas from ~kpc scales down to the BH is indeed ef- 
ficient, as discussed in §2. 1-2.2. Feedback is modeled 
not as a full description of the outflow, but by transfer- 
ring a fraction of the radiative output of the BH to the 
thermal energy of the surrounding gas. The feedback 
energy output is thus 



"eed 



(11) 



where e,. ~ 0. 1 is the standard radiative efficiency, and 
6/ is the fraction of the radiative output that couples 
to the gas. Such models naturally yield a relationship 
between Mbh and Mjph, and can reproduce the evolu- 
tion with redshift as observed from broad-line quasar 
samples (e.g.. Booth & Schaye 201 1). To reproduce 



the observed normalization for these relations requires 
6/ ^ 0.05-0.15 (depending on the specific subgrid pre- 



scriptions for feedback; Springel et al. 2005b ; Di Matte"o| 
[eFal. 2008 ; Booth & Schaye 2009; Teyssier et al. loTT) 
suggesting relatively efficient coupling between the BH 
luminosity and the surrounding gas. A detailed study of 



the output of one model ( .Booth & Schaye 2010 1 sug- 
gests that the most fundamental relationship is between 
Mbh and the binding energy of the host dark-matter 
halo. In this model, the BH grows until it is massive 
enough to unbind gas not only from its host galaxy but 
from the entire surrounding halo, and it is this feedback 
(along with correlations between the galaxy and Mhaio) 
that give rise to the observed BH-spheroid correlations. 
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Recent observational results suggest that such a Mbh- 
Mhaio relation is not fundamental for lower-mass AGN 
(e^ lGreene et al.| [SOOS] [20T0l [Kormendy & Bender] 



201 1\ but it remains a valid possibility for higher-mass 



objects Mbh ?J 10 Mq for which BH, galaxy spheroid, 
and dark-matter halo correlations are more difficult to 

Son) . 



disentangle (jVolonteri et al. 



While feedback has naturally occupied a great deal of 
attention in understanding BH-spheroid relationships, 
recent work has also proposed that these relations do not 
imply a causal link between BHs and their host galaxies 
or halos, but are simply a consequence of the hierar- 
chical growth of structure (e.g., |Peng[ |2007[ |Jahnke &| 
Maccio 201 IV If mergers of galaxies are accompanied 



by the corresponding mergers of their central BHs, then 
a succession of mergers will tend to drive the result- 
ing galaxies toward the "average" relationship between 
BH and spheroid mass, even if the initial galaxies show 
little or no correlation between BH and spheroid prop- 
erties. These processes may also help explain the in- 
creased scatter at the low-mass end of the BH-spheroid 
relationship, for which the galaxies may have experi- 
enced fewer mergers. While in this picture mergers play 
the predominant role in producing the BH-spheroid re- 
lationships, AGN feedback may be important in setting 
the relative rates of BH and galaxy growth over the his- 
tory of the Universe. 

5.2.3. The impact of feedback on gas in galaxies and 
halos 

The energy released by growing BHs can impact not 
only on the stellar content of their host galaxies (as sug- 
gested for feedback models of the BH-spheroid relation- 
ships) but also the diffuse gas in their host dark matter 
halos. In massive (Mhaio lO'"^ Mq) halos with hot, 
ionized atmospheres, heating from relativistic jets (as 
discussed above) can regulate the temperature of the ha- 
los and prevent cooling to form new stars, as required in 
order to produce the observed stellar mass and luminos- 
ity functions of galaxies (e.g.. Bower et aT) |2006[ |Cro-| 



ton et al.||2006[|Bower et al.||200 8). Mechanically dom- 
inated AGN activity is observed to be most prevalent in 



passive systems in such massive halos (e.g.. Best et al. 



Hickox et al.' '2009'; Wak e et al.| |20081 [MandeH 
baum et al. , 2009, Smolcic, 2009] ); see §3.4 & Fig. 5. In 



contrast, moderate-luminosity radio-quiet AGN appear 
to have very little impact on the stellar populations or 



star formation rates of their host galaxies (e.g., Hickox 



leFaLl [2009] [Cardamone et"aLl [20T0l [Xue et al.| |2010r 
ILutz et al.' '20T0| |Shao et al.| [20T0l [Ammons etal] 
[201 1 , Mullaney et aH|201 lb |), although outflows in such 
systems may still play a significant role in moderating 



the growth of the BH (e.g., |King|2010b| Mullaney et al., 
in prep). 

Some of the most compelling evidence for AGN feed- 
back comes from observations of the hot gas itself, 
which represents the majority of the baryonic content 
of group- and cluster-scale systems. Of particular in- 
terest is the luminosity-temperature {Lx-Tx) relation, 
which is broadly related to the baryon fraction in a halo 
as a function of Mhaio. Observations show that the Lx- 
Tx relation is remarkably steep at the low-temperature 
end, suggesting that small to moderate-mass groups 
have lower baryon fractions in the centers of their dark- 
matter halos than higher-mass systems (e.g., |Sun et al. 



[20091 [Giodini et aT) [2009J. Outflows must have re- 
moved some of the low-entropy gas from the centers 
of these halos, and simulations indicate that winds from 
stellar processes (stellar winds and supernovae) are in- 
suflicient to produce this deficit. However, models sug- 
gest that feedback from growing massive BHs (indeed, 
the same outflows that are invoked to produce the BH- 
spheroid relation) can naturally produce the Lx-Tx rela- 



tions similar to those observed (e.g.. Bower et al. 



Boofli & Schaye 2010 McCarthy et al. 2010 



J2008 
20TT] i. 



Unlike models for the galaxy luminosity function or 
BH-spheroid relationships, these predictions are inde- 
pendent of assumptions about the star formation pro- 
cesses within the host galaxy and depend only on the 
coupling of the AGN to the surrounding gas. As such, 
they represent one of the most promising avenues for 
quantifying the impact of growing BHs on their host 
galaxies and halos. 



6. Concluding remarks 

In this review we have discussed the processes and 
mechanics of BH growth and explored the connections 
between the growth of BHs and star formation, host- 
galaxy properties, signatures of interactions/mergers, 
and environment. We have explored these processes in 
a very broad range of systems, from the closest AGNs 
to AGNs at the highest redshifts, and also the formation 
of seed BHs. Our overarching aim has been to address 
the question: What drives the growth of black holes ? 

Surprisingly, there are often few external morpholog- 
ical and kinematical clues that reveal when significant 
BH growth has been triggered, with AGN host galaxies 
showing similar overall properties to inactive galaxies 
of comparable stellar mass. This appears to be true even 
down to small ^ 10-100 pc scales. However, the accu- 
rate identification of gas inflow on these scales is mostly 
restricted to low-luminosity AGNs in the local Universe, 
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and significant differences may be present in the cen- 
tral X! 10-100 pc regions of luminous AGNs. Compa- 
rable linear spatial scales should be achievable for dis- 
tant and high-luminosity AGNs over the coming years 
with the final configuration of ALMA at submm-mm 
wavelengths, which will be provide direct insight into 
the properties and extent of the molecular gas and dust- 
continuum emission on ^ 100 pc scales. A clear compli- 
cation in identifying any large-scale catalysts that may 
have triggered AGN activity is the timescale required 
for the gas to flow from the host galaxy down to the 
BH, which can be » 10^ yrs. Over such long timescales 
any clear morphological or kinematical signatures of the 
triggering event may be lost and therefore care needs to 
be taken when assessing the role of galaxy major merg- 
ers or interactions in the fueling of AGN activity. 

A clear dependence between the SFR and luminous 
AGN fraction is found in both distant and local galaxies, 
implying a connection between BH and galaxy growth. 
This finding is perhaps not unsurprising since both star 
formation and AGN activity are predominantly driven 
by a cold-gas supply, either from the host galaxy and 
larger-scale environment, or from stellar winds and cir- 
cumnuclear supernovae events. However, these results 
may also indicate a more direct and less coincidental 
connection between star formation and AGN activity in 



at least some sytems (such as in the common evolution- 
ary scenario shown in Fig. 6). High-resolution obser- 
vations with ALMA and JWST of well-selected galaxy 
samples could directly test this evolutionary scenario by 
simultaneously providing SFRs, gas-inflow rates, and 
BH mass accretion rates as a function of stellar age. 

Recent results have shown that environment plays a 
major role in the triggering of AGN activity. The global 
peak in AGN activity in overdense regions occured at 
higher redshifts than in the field, with a corresponding 
more rapid decline to the present day. The most rapid 
BH growth also appears to be found in dark-matter ha- 
los of « lO'^-lO'-' Mq (see Fig. 5), with more massive 
dark-matter halos found to harbour low-accretion rate 
radio-loud AGNs hosted by relatively quiescent galax- 
ies. This mass dependence on AGN activity may be due 
to "quenching" of the cold gas in massive dark-matter 
halos, which could be largely assisted by large-scale 
mechanically powerful (but often radiatively weak; see 
Fig. 7) outflows, such as those seen in nearby galaxy 
clusters. The strong dependence on the fraction of 
galaxies that host radio-loud AGN with galaxy mass 
provides further support for this picture. 

However, on the basis of the current results it is not 
immediately apparent why there exists a tight correla- 
tion between the mass of the BH and the mass of the 
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spheroid in nearby galaxies. Many models argue that 

energetic radiatively driven outflows (such as those il- 
lustrated in Fig. 7) regulate star formation in the host 
galaxy and "forge" the BH-spheroid mass relationship. 
Energetic radiatively driven outflows are clearly seen in 
many luminous AGNs and may be a ubiquitious feature 
of the AGN population but it is not yet clear what im- 
pact these outflows have on star formation in the host 
galaxy. Spatially resolved spectroscopy of AGNs with 
IFU observations may provide important insight on this 
issue since they can reveal the presence of large-scale 
outflowing gas and address when and where energetic 
outflows are found. However, given the relatively mod- 
est BH mass accretion and SFRs for the majority of the 
AGN population (e.g., see Fig. 4 for the AGN fraction 
as a function of different AGN selection techniques), 
a "regulatory" mechanism may not be necessary every 
time the BH and galaxy grows. It is presumably only the 
rapidly evolving and gas-rich systems where the BH and 
galaxy growth can be so high that they can move sub- 
stantially around the BH-spheroid mass plane on short 
timescales. These rapdily evolving systems are there- 
fore the most promising sites to search for the evidence 
of large-scale energetic outflows. We should also bear 
in mind that the tight BH-spheroid mass relationship is 
only apphcable for galaxies hosting classical bulges - 
there are no significant BH-galaxy correlations for sys- 
tems hosted in pseudo bulges. Since galaxy major merg- 
ers are likely to be the catalyst for the formation of clas- 
sical spheroids, this suggests that galaxy major mergers 
played a major role in forging the BH-spheroid mass 
relationship. 

Although rare, luminous AGNs such as the quasars 
selected in optical surveys may play a much more im- 
portant role in the evolution and growth of BHs and 
galaxies than would be apparent from their relatively 
low space density. An example of this is the finding 
that the most luminous AGNs in the local Universe are 
responsible for « 50% of the integrated BH growth, de- 
spite only comprising » 0.2% of the optical AGN popu- 
lation. However, since AGNs can be extremely variable, 
any measurement of the BH mass-accretion rate only 
provides a snapshot of its growth path (e.g., see Fig. 5 
of Hopkins & Quataert 2010, which shows factors of 
« 100 variation in the predicted gas inflow rate over Myr 
timescales). This adds significant complications to the 
interpretation of individual objects, which can be over- 
come to some extent from studying AGN populations as 
a function of a more stable property than mass accretion 
rate such as BH mass, steUar mass, and star-formation 
rate. This should be possible over the coming years with 
the advent of sensitive large facilities, particularly those 



that combine depth with breadth and yield large source 

densities (e.g., BigBOSS, LSST, and e-ROSITA). 

The current theoretical models focus on either small- 
scale processes at high resolution, with Umited con- 
straints on larger-scale processes, or explore the large- 
scale growth of BHs and galaxies in cosmological simu- 
lations using sub-grid models to account for the smaller- 
scale physics. Both approaches have clear strengths 
but evident weaknesses (i.e., either limited spatial and 
time resolution or limited cosmological volume). Since 
AGNs are gas-driven systems fed from gas inflow down 
to < 10 pc scales from typically > 1 kpc scales, high- 
resolution hydrodynamical simulations over a broad 
range of size scales would be particularly beneficial for 
a detailed understanding of AGN fueling and feedback 
processes. More detailed theoretical insight may also 
come from improvements in the sub-grid models using 
the physics learnt from yet higher-resolution small-scale 
models to produce more reahstic fuehng and feedback 
prescriptions for the cosmological simulations. 

The next 5-10 years therefore promise huge ad- 
vances in our understanding of "What drives the growth 
of black holes?". We look forward to hosting another 
exciting workshop and reviewing our progress on this 
issue over the next decade! 
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